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Abstract 

The linearized dynamic and kinematic equations for the 
three-dimensional motions of a faired towing cable in a 
uniform flow are formulated with special attention to the 
hydrodynamic forces. Sinusoidal kinematic boundary conditions 
are applied to the cable at the tow point and simple dynamic 
boundary conditions are applied at the towed body. Analysis 
is done in the frequency domain. The resulting equations are 
solved numerically. The method is shown to be applicable over 
a range of practical cases. A computational difficulty 
involving the interaction of growing solutions with round-off 
error is identified. Demonstrably the difficulty may be put 
off by increasingly precise computation. The effect of import¬ 
ant physical parameters of the problem on the motions trans¬ 
mitted from the tow point down the cable to the towed body are 
evaluated. Results are understood in terms of the natural 
frequencies of the oscillatory modes which combine to produce 
the cable dynamics. A design principle is established for 
minimizing towed body motions by partially decoupling the tow 
point from the towed body using moderately extensible cables. 


Thesis Supervisor; J. Nicholas Newman 
Title: Professor of Naval Architecture 


2 












ACKNOWLEDGEMENTS 


I greatly appreciate the advice of my thesis supervisor. 
Professor J. Nicholas Newman and the efforts of the other 
members of my committee. Professors E. E. Covert, J. H. Milgram 
and J. T. Christian. Financial support for this work was 
provided by the Naval Resear c h and Development Center, 
Carderock, Maryland,through the ef£orts“ofTlr. W. E. Smith, 
Head , Shi p Dyn^ics Simulatipn^Branch. Captain W. F. Searle, 
TJSNCret) was influential in the selection of the topic and in 
advising me concerning my doctoral studies and deserves 
special thanks. The capable typing of Miss Sandra Vivolo was 
a considerable assist. 

The value of the patience and understanding of my wife, 
Barbara, through this whole experience is inestimable. 


3 






I 





I 

i 






TABLE OF CONTENTS 

Title Page 1 

Abstract 2 

Acknowledgements 3 

Table of Contents 4 

List of Figures 6 

List of Tables 7 

List of Symbols 8 

1. Introduction 14 

Background 14 

Antecendents 16 

2. Formulation 20 

Coordinate Systems 20 

Dynamic Equations 22 

Kinematic Equations 28 

Lagrangian Coordinate 30 

Non-dimensionalization 33 

Boundary Conditions 35 

3. Solution 41 

Frequency Domain Analysis 41 

Steady Configuration Solution 47 

Dynamics Solution 50 

Solution Difficulties 54 

4. Results 58 

Basic Cable System 60 

Variation of Cable Parameters 62 

Natural Frequencies 71 

4 




1 


i 







5. Conclusions 


77 


Design Impact 

Suggestions for Future Work 
Appendices 

I. Hydrodynamic of Faired Cables 
Introduction 
Added Mass 
Loading Functions 
Thin Wing Theory 
Results 

II. Hydrodynamics of Towed Body 

III. Special Cases 

Neutrally Buoyant Cable Without Drag 

Neutrally Buoyant Cable and Body With 
Normal Drag 

IV. Trial Cable 

Bibliography 

Biographical Note 


77 

78 

79 
79 
82 
86 
91 

100 

102 

106 

108 

111 

115 

116 
121 


5 




# 





LIST OF FIGURES 


2.1 Coordinate Systems 23 

2.2 Differential Element of Faired Cable 23 

3.1 Typical Steady Configuration 

4.1 Motions and Tensions—Basic Cable 61 

4.2 Vertical and Horizontal Motion—Length Variation 64 

4.3 Lateral Motion and Tension—Length Variation 65 

4.4 Horizontal Motions—Steady Velocity Variation 66 

4.5 Vertical Motions—Steady Velocity Variation 67 

4.6 Lateral Motions--Steady Velocity Variation 68 

4.7 Horizontal and Vertical Motions—Elasticity Variation ^6 

4.8 Horizontal and Vertical Motion with Cj^ = 0.0 ^2 

4.9 Lateral Motions and Tension with Cj^ =0.0 

4.10 Horizontal and Vertical Motions—Very Elastic, Cj^=0.0 75 

I-l Typical Fairing Segment 83 

1-2 Rotational and Translational Velocities of Fairing 

Segment 83 

1-3 Inflow Velocity and Components 89 

1-4 Lateral Velocity and Force Diagram 89 

1-5 Fairing Dimensions 89 


6 



I 



LIST OF TABLES 


3.1 

Comparison of Single and 
Programs 

Double Precision 

57 

4.1 

Changes from Basic Cable 
Length Variations 

Parameters 

for 

63 

4.2 

Changes from Basic Cable 
Tow Velocity Variations 

Parameters 

for 

63 

4.3 

Changes from Basic Cable 
Elasticity Variations 

Parameters 

for 

63 


7 




LIST OF SYMBOLS 


cross-section area of unloaded cable 

cross-sectional area of towed body 

non-dimensional body displacement defined 
on pg. 33 

accelerations in the x, y and z directions 

fairing, hydrodynamic loading function, sine 
series coefficients 

fairing breadth in the loaded condition 

effective fairing breadth in unloaded condition 

fairing, hydrodynamic loading function^ cosine 
series coefficients 

non-dimensional body drag defined on pg. 39 

fairing drag coefficient 

drag coefficient of towed body 

fairing chord length in the loaded condition 

effective fairing chord length in unloaded 
condition 

non-dimensional drag defined on pg. 33 
non-dimensional drags defined on pg. 34 

7 

non-dimensional cable and fairing displacement 
per unit length defined on pg. 33 

modulus of elasticity of the cable 

external force on fairing 

amplitude of external force on fairing 

normal and tangential cable loading functions 

towed body disturbing forces 

hydrodynamic function defined on pg. 95 




acceleration due to gravity 

hydrodynamic function defined on pg. 95 

non-dimensional fairing mass per unit length 
defined on pg. 33 

heave perturbation velocity in x , y , z 
coordinate system o o o 

mass moment of inertia of fairing per unit 
length about y axis 

Bessel's functions of the first kind defined 
on pg. 98 

non-dimensional towed body added mass defined 
on pg. 39 

length of the cable 

non-dimensional towed body mass defined on 
pg. 39 

mass of towed body 

hydrodynamic moment derivatives per unit 
length for moments about y axis with respect 
to subscripted variable 

mass per unit length of cable and fairing 
in loaded condition 

mass per unit length of fairing and cable in 
unloaded condition 

fairing mass per unit length in loaded 
condition 

effective fairing mass per unit length in 
unloaded condition 

non-dimensional steady tension defined on pg. 33 

non-dimensional perturbation tension defined 
on pg. 33 

complex amplitude of n 

real and imaginary parts of complex 
amplitude of n 

magnitude of n 

9 





Pf q/ r 

rotational velocities about the x, y and z 
axes respectively 

s 

Eulerian distance coordinate along the cable 

6U 

surge perturbation velocity in x^, y^, 
coordinate system 


sway perturbation velocity in x , y , z 
coordinate system 

T 

tension in the cable 

T* 

steady configuration tension in the cable 

/S 

T 

perturbation tension in the cable 

t 

time 

U 

o 

steady tow velocity 

u, w 

velocities in the x, y, and z directions 
pg. 14 to 32 

U, V, w 

non-dimensional perturbation velocities in 

X, y, and z directions pg. 33 to 115 

U*^ V*, w* 

steady configuration velocities in the 

X, y, and z directions 

/S ^ ^ 

U, V, w 

perturbation velocities in the x, y, and z 
directions 

U, V, w 

complex amplitudes of u, v, w pg. 33 to 115 

’^R, ’^R' 

real and imaginary parts of complex 
amplitudes of u, v, w 

H 

H 

> 


U , V , w 

o o o 

velocities in the x , y , z coordinate 
directions 


u V , w magnitude of u , v , w 
om om om ^ o o o 


V 

volume per unit length of fairing and cable 
in loaded condition 

V 

volume per unit length of fairing and cable 
in unloaded condition 


volume of towed body 

X, Y, Z 

force components in the x,y, and z directions 
with superscripts indicating contributions 
from various effects 

10 





I 



hydrodynamic force and moment derivatives per 
unit length of cable fairing in the x and z 
N* directions and about the x, y, and z axes 
^ respectively 

non-dimensional body disturbing forces defined 
on pg. 39 from pg. 39 to pg. 78 

cable oriented Cartesian coordinate axes 

locally vertical Cartesian coordinate axes 

distance from y axis to center of gravity of 
fairing segment 

effective distance from origin of cable 
coordinate system to center of gravity of 
fairing segment 

distance from y axis to point through which 
resultant viscous drag force of fairing per 
unit length acts 

Bessel's functions of the second kind defined 
on pg. 98 

composite z component of fairing gravitation and 
inertial forces combined with external force of 
cable on fairing 

hydrodynamic force derivatives per unit length 
for forces in z direction with respect to 
subscripted variable 

z component of viscous drag force per unit 
length 

z component of force on fairing from cable 

amplitude of fairing motion in z direction 
pg. 79 to pg. 101 

fairing rotation angle about the y axis 
pg. 79 to pg. 101 

amplitude of fairing rotation angle pg. 79 
to pg. 101 

non-dimensional surge, heave and sway velocities 
in the cable coordinate system pg. 35 to pg. 73 

complex amplitudes of a, B, e pg. 35 to pg. 73 



















































angle of inclination of inflow velocity vector 
pg- 88 to pg. 89 

non-dimensional modulus of elasticity of the 
cable defined on pg. 33 

non-dimensionalization factor defined on pg. 33 

non-dimensional Lagrangian distance along the 
cable defined on pg. 33 

non-dimensional lateral added mass defined on 

pg* 33 

reduced frequency defined on pg. 95 

non-dimensional transverse added mass defined 
on pg. 33 

added mass of towed body 

non-dimensional angular frequency of excitation 

non-dimensional distance of fairing center of 
gravity from cable coordinate system origin 
defined on pg. 33 

water density 

Lagrangian distance coordinate along the cable 

non-dimensional time defined on pg. 33 

steady configuration trail angle in non- 
dimensional equations 

critical trail angle 

trail angle pg. 14 to pg. 32 

perturbation trail angle pg. 33 to pg. 115 

steady configuration trail angle 

perturbation trail angle 

complex amplitude of ^ pg. 33 to pg. 115 

real and imaginary part of complex amplitude of cj) 

kite angle pg. 14 to pg. 32 

perturbation kite angle pg. 33 to pg. 115 



Ip* 

$ 

¥ 

h' h 


(0 

0 ) 

0) 

w 

w 

U) 


Pt 


P5, 


Tt 


TJl 


steady configuration kite angle 
perturbation kite angle 
conplex amplitude of ip 

real and imaginary parts of complex amplitude 
of ip 

excitation frequency in radians per sec. 
transverse pendulum natural frequency 
lateral pendulum natural frequency 
transverse vibrating string natural frequency 
lateral vibrating string natural frequency 
longitudinal stress vibrations natural frequency 


Other symbols are used only locally and defined where used. 


13 


I 

I 


I 



chapter 1 
INTRODUCTION 


Background 

A ship tows a body suspended at the end of a long cable. 
The cable is slender and curvatures are small so bending 
moments are negligible. The ship proceeds at velocity in 
a seaway and bounces around a lot. As it moves along its 
motion is somehow transmitted to the towed body through the 
cable. The central question addressed by this thesis is, 
"what motions are transmitted to the towed body through the 
cable?" Real towed body systems now generally employ a 
fairing that is mostly free to rotate about the cable. So 
the cable is assumed to have a freely rotating fairing 
attached. Real cables stretch under load so the cable is 
assumed to be extensible. 

In the absence of any excitation a steady towing situa¬ 
tion exists, a problem solved by Pode ( 1951 ). The basis of 
the solution of the unsteady problem with an excitation 
present will be the assumption that the motion consists of 
small perturbations of the motion variables about their 
steady values leading to linear equations. These equations 
may be solved in the time domain with the method of charact¬ 
eristics or in the frequency domain assuming sinusoidal 
excitations of the tow point. Although both approaches have 
merit the frequency domain solutions are compatible with the 
existing body of knowledge concerning ship motions and with 
a formulation of the hydrodynamics of cable fairings derived 
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from thin wing theory. Knowledge of the ocean wave spectrum 
combined with response amplitude operators for a given tow 
ship can be converted to a spectrum of tow point velocities 
making an understanding of cable dynamics in the frequency 
domain appropriate. The frequency domain approach will be 
pursued here. 

Before proceeding to the details of this paper it may be 
instructive to compare the problem posed to the classical 
problem of the motions of a vibrating string. A string is 
stretched between two points whose locations are fixed. The 
tension in the string is considered constant throughout its 
length. For certain conditions its rest position is a 
straight line between the end points. The string is considered 
to be inextensible. Its motions are described by the wave 
equation \p = c^ijj where ip is the displacement of the string 
from its rest position. 

The problem posed in this thesis differs from the class¬ 
ical vibrating string problem in the following respects: the 
boundary conditions are more complicated than the simple kine¬ 
matic conditions described above. One boundary of the tow 
cable is subject to a kinematic condition consisting of pre¬ 
scribed motions. The other boundary is subject to a dynamic 
condition consisting of prescribed forces, themselves functions 
of the kinematic variables at the boundary. The tension in the 
cable is not constant throughout. In fact the tension increases 
in some fashion from the body end to the towpoint end of the 
cable. The equilibrium condition of the cable is a non-linear 
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space curve. The cable is extensible, and finally it is 
operating subject to hydrodynamic and gravitational forces. 

The towing cable is a complicated combination of a vibrating 
string, a compound pendulum and a distributed- mass-spring 
system all with significant viscous damping. 

The solution of the cable problem will in various ways 
include the solutions to these classical problems. Of course, 
the coupling between the various problems is considerable and 
so the features of each are only outstanding in certain 
special cases. This thesis will formulate the problem in its 
complexity, apply the frequency domain analysis to plausible 
designs and bring the above notions to mind again while 
interpreting the results. 

Antecedents 

The problem of determining the cable configuration for a 
steady tow was comprehensively solved by Pode ( 1951 ). Both 
before and after Pode's solution investigators have studied 
hydrodynamic loading functions for cables. Their work is 
summerized by Casarella and Parsons ( 1970 ). In the same 
survey a number of investigations into the unsteady motions 
of cables are also summarized, of which the work of Whicker 
is very relevant to this thesis. 
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Whicker ( 1957 ) studied the dynamics of an unfaired, 
inextensible cable with drag but no added mass. He examined 
the transverse and longitudinal motions of the cable in a 
plane parallel to the towing direction. He presented two 
schemes for solution, one employing the method of character¬ 
istics and the other a linearized frequency domain approach. 
Neither solution scheme could be fully implemented without 
extensive computer aid so only simplified special cases were 
worked out. The special cases considered in this thesis in 
Appendix III where coincident with those of Whicker, verify 
his results. 

Abkowitz ( 1965 , 1967 ) considered the torsional and 
lateral vibrations of a towing cable with the fairing bound 
to the cable. His interest was in establishing cable para¬ 
meters for which this system would be stable with respect to 
hydrodynamic flutter. Abkowitz showed that for some cable 
and fairing parameters in the range of practical interest 
stability was marginal and flutter could be expected. Pri¬ 
marily because of flutter and more general kiting instabilities 
most faired cable systems now employ fairings which are more 
or less rotationally unrestrained. 

Schram and Reyle ( 1968 ) used the method of characteristics 
to analyze the transverse, longitudinal and lateral motions 
of an inextensible cable system with drag but no added mass. 
They developed transfer functions comparing the amplitude of 
surface excitations to the amplitude of towed vehicle motion. 
Their analysis was complicated by the assumption that the 
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cable was inextensible. Then a pair of equations in the first 
order system are parabolic instead of hyperbolic creating 
difficulties for the well-posedness of the initial value 
problem. 

Huffman (1969) considered an extensible bare cable 
towed system undergoing transverse and longitudinal vibrations 
while being towed through air. He used the method of char¬ 
acteristics to describe cable configurations, variations in 
tension and towed body displacements from equilibrium for 
various combinations of cable parameters. His method was 
fully implemented numerically and since all the equations 
were hyperbolic, was well posed. As mentioned before however, 
these solutions in the time domain are not particularly 
compatible with the spectral analysis of ship motions. Of 
course, Huffman's interest was not in sea systems but in air 
towed systems so that this incompatibility was not a detriment 
to his selection of the method of characteristics as the 
proper approach to the problem solution. 

Kerney (1971) used a small pertubation scheme to analyze 
the two-dimensional motions of an inextensible cable with 
added mass. Kerney's frequency domain analysis differs little 
from the method applied here or from the original proposal of 
Whicker. In implementing his solution numerically he encoun¬ 
tered difficulties left unresolved. 

The concern of the present study is an extensible cable, 
with a freely rotating fairing attached, towing a body and 
excited by transverse, longitudinal and lateral motions 
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emanating from the tow ship. The significance of the freely 
rotating fairing is that it greatly reduces drag on the cable. 
It can transmit no moment to the cable but it can transmit a 
side force the nature of which will be investigated. The 
results of the study concern the response of the cable and 
the towed body to sinusoidal excitations. 
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Chapter 2 
FORMULATION 


Coordinate Systems 

A right-handed Cartesian coordinate system moves with the 
local position of the cable with its axis along the 
direction of the mean tow velocity U^. The y^ axis points 
opposite to the direction of gravitational attraction and the 
axis points to starboard. Another right-handed Cartesian 
system unsubscripted is imbedded in the towing cable in such 
a way that were the tow ship stationary with the cable hanging 
straight below this x, y, z system would point in the same 
directions as the x^, y^, z^ system at all points along the 
cable. When the tow ship is moving with velocity through 
the stationary ocean the cable streams behind and below the 
tow ship in a complicated space curve so that at each point 
along the cable the cable coordinate system will point in 
different directions from the x ,y , z system. The relation 
between the x, y, z system and the x^, y^, z^ system is 
defined at each point of the cable by two angles, the trail 
angle (}) and the kite angle ip. The trail angle (}) is the angle 
measured between the y axis of the cable system and the x^ 
axis of the ship system considering a rotation about the z 
axis. The kite angle ip is then formed by a rotation about 
the X axis. The relationship between the two systems shown 
in figure 2.1 is given by the following transformation matrix 
which is formed by a procedure similar to that of Abkowitz 
(1969) . 
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( 2 . 1 ) 


X 


sin(j) -cos4) 0 


X 

o 

y 

= 

cosi|< cos(i) cosi|^ sin4) sini|^ 


^o 

z 


-sini|) cos(l> -sini|^ sini}) cosip^ 


■^o> 


In the case of steady towing as shown by Pode (1951) 
there is no kiting but there is a trail angle (p. In the 
ensuing analysis perturbations about the steady values are to 
be used. The trail angle and kite angle are written as the 
sum of steady values and perturbation values which are small 
compared with the steady values, thus 

({) = (()* + $ 

ip = ip* + ip and furthermore 4^* = 0. 


Under these circumstances the transformation matrix elements 
may be expanded in a Taylor Series giving for the zeroth and 
first order in the perturbation quantities 


( > 
X 


sin4>*+cos(|)*$ -cos4>*+sin(|)*$ 0 


X 

O 

y 

= 

cos4)*-sin<})*$ sin4)*+cos<)*$ ip 


^o 

lZj 


cos<p* -ip sinp* Ij 

' 



( 2 . 2 ) 
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Angular velocities to first order in perturbation quantities 
are 9il)/9t and 9$/9t. Accelerations in the x, y, and z 
directions are given by 


a 


X 


9t 





a 


z 


9t 


V* 


9t 


{2.3a-c) 


where u, v, and w are the velocities of the cable in the x, 
y, and z directions respectively. As before the velocities 
are the sum of steady values and perturbation values, and can 
be written as 


u = u* + (i 


V = V* + V 

w = w* + ^ but w* = 0 so w = w 

Dynamic Equations 

Consider a differential element of a faired towing cable 
such as shown in figure 2.2. The components of forces in 
each of the three directions of the cable coordinate system 
must be in equilibrium. For the unsteady case contributions 
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Figure 2.1 
Coordinate Systems 



Figure 2.2 

Differential Element of Faired Cable 
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to the forces are made by the inertia of the cable and the 
fairing, the tension in the cable, the weight of the cable 
and the fairing and the fluid flow over the surface of the 
fairing. Assume that the fairing is segmented into very 
short segments each of which is mounted on frictionless 
bearings on the cable so that forces are transmitted from the 
fairing to the cable but no twisting moment is transmittf^d. 
This assumption appears to be a fairly accurate description 
of many real systems. 

The inertial forces acting on a differential element are 


X - -mds - V* 



(2.4a-c) 



where the force components in the respective directions are 
given by the upper case letters. The mass per unit length 
of the cable and fairing is m and the length of the differ¬ 
ential element in the loaded condition is ds. 

The gravitational forces are: 



= -(m-pV)ds g (sin(|)*+cos^*$) 


(2.5a-c) 


Z = (m-pV)ds gipsin^* 
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where p is the density of the water in which the system is 
operating, V is the volxime per unit length of the fairing and 
cable in the loaded condition and g is the acceleration due to 
gravity. 

The tension forces acting on the cable are 
x(T) = .T. II* as - i II* ds - T* II ds 


(T) 

9s 


ds + ^ ds 


(2.6a-c) 


Z™ = T* ||dS 

I d S 


where T is the tension which is the sum of a steady value and 
a perturbation value, that is T = T* + T. 

The hydrodynamic forces acting on the cable may be 
written as 


X 


(h) ^ j^*(h) ^ JJ(h) 


y(h) ^ y,(h) ^ y(h) 

z(h) = s.*(h) + z(h) 


where as before, since the steady motion of the cable is in 

* (h) 

the x-y plane, Z = 0. Appendix I contains the details of 

the hydrodynamic analysis and shows that the following 
equations may be written for (2.7a-c). 
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= “Cpip cds (a^+a 2 Sin(})*) sin(t)*U^ 

-Cj^'^p cds [ { (a^+a 2 sin 4 )*) 2sin^4)*+ (aj^+ 2 a 2 sin(})*) 

*\ /V 

cos^4)*}U^u+a, sin(j)*cos(j)*U <r]-^p7rb^ds (2.8a) 

OX O d t 

= -Cj^^p cds (bj^+b 2 COS(j)*) cos<t»*U^ 

-Cjjip cds [b^sin(})*cos4)*U^u+{2 (bj^+b 2 COS(})*) cos^^* 

+ (b^+ 2 b 2 COS(j)*) sin^(J)* }U^v] -^pirb^dsu*^ (2.8b) 


= “ipCjj cds (a^+a 2 sin(()*)U^w - ^irpc^ds 


-4m 


V* 3i 


ds 


8t 


(2.8c) 


where b, c, m, x_, C^, a,, a_, b, and b„ are defined in 

Qj X) X z 1 Z 

Appendix I. Equation (2.8c) is restricted to motions of the 

cable whose frequency corresponds to low reduced frequency, 

coc 

— and for fairings made of materials with densities close to 
^o 

that of sea water. 

Combining equations (2.4a-c), (2.5a-c), (2.6a-c), and 

(2.8a-c) gives the dynamic equations for a differential 
element of cable. 

-mds (•l^v*'|4) + (m-pV) ds g (cos4)*-sin4)*$) -T* |X*ds-T|^*ds-T*|i ds 

dtldt dS dS dS 

-Cj^^pcds (a^+a 2 Sin 4 »*) sin4)*U^-Cjj^pcds [{ (a^+a2sin4)*) 2sin^(j)* 

+ (aj^+2a2sin<j)*) cos^4)* }U^(i+a^sin<t)*cos4i*U^'(^]-iJ-pTrb^ds = 0 (2.9a) 
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(bj^+b 2 COS(j)*) cos4)*U^-Cplpcds [bj^sin<J)*cos(j)*U^u+{2 (bj^+b 2 Cos(j)*) 



(2.9b) 



-^pCpCds (a^+a2sin4)*) U^w = 0 


(2.9c) 


Equations (2.9a,b) contain terms which are steady and terms 
which are first order in the time varying perturbation 
quantities. The steady terms may be separated from the time 
varying terms to yield a pair of equations describing the 
steady motion of the cable and a pair which pertain only to 
the time varying parts of the motion. The steady equations 
are 

(m-pV)ds g cos4)*-T*|^ ds-ipC^cds (aj^+a 2 sin(|)*) sin(|)*U^=0 


(2.10a) 



(2.10b) 


The time-varying equations are 


-mds (•^ -V* 


|^)-(m-pV)ds g sin<J)*$-T ds-T* ds 

o t d S d S 


-Cp^pcds [ { (aj^+a2sin(|)*) 2sin^4)*+ (aj^+2a2sin4)*) cos^4)* }U^u 


+a^sin(|)*cos(t)*U^<^] “i^piTb^ds = 0 


(2.11a) 
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-mds ("l^ +u*|^)-(m-pV)ds g cos<|)*$+|^s-Cj^ipcds [b^sin(})*cos<j)* 
U^u+{2 (bj^+b 2 COS(t)*) cos^<f)*+ (bj^+2b2) cos4)*) sin^cf)* }u^v] 


-^PTTb^ds u* = 0 


(2.11b) 


Equations (2.10a,b) are the classical steady equations for a 
towed cable as determined by Pode (1951) and others. Equations 
(2.11a,b) and (2.9c) are the dynamic equations for the 
linearized three-dimensional motions of a faired cable. 

Kinematic Equations 

Kinematic relations for the differential element of 
cable must also be satisfied. Again consider the element 
shown in figure 2.2. Points P and Q can move relative to 
each other in three ways, one way along each of the intrinsic 
coordinate directions. Motion along the y direction comes 
from stretching the element. Motion in the x and z directions 
is caused by relative rotation of the element. The following 
relations express these relative motions. In the y direction 



(2.12a) 


In the X direction 



(2.12b) 
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and in the z direction 



(2.12c) 


As in the case of the dynamic equations, the relations of 
equations (2.12a,b) contain steady terms and time-varying 
terms. The equations may then be separated each yielding two 
equations, one for the steady motion and one for the time- 
varying motion. Since there is no relative motion between 
points P and Q for the steady motion of the cable, the trivial 
results of the steady portion of equations (2.12a,b) is 
expected, namely that they are satisfied by u* = sincf)* and 
V* = cos(p*. The time-varying portions give 


ds + u* ds + u ^ ds = (ds) 

oS oS oS dt 


(2.13a) 



(2.13b) 


The dynamic equations and the kinematic equations for the 
time-varying three-dimensional motions of the cable are six 
equations in six unknowns, namely the three perturbation trans¬ 
lational velocities u, v, w, the two perturbation angular 


velocities $ and ip and the perturbation tension T. The 
equations for the x and y motion are coupled together and de¬ 
coupled from the side motions in the z direction. All the equa¬ 
tions are written in the two independent variables s and t. The 
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equations have space varying coefficients in the variables 
u*, V*, <l>* and T* a set of variables which trivially reduce 
to (})* and T*. These quantities may be determined from the 
solution of the steady equations (2.10 a,b). Before pro¬ 
ceeding to methods of solution some additional operations may 
be performed on the cable equations to simplify the analysis. 


Lagrangian Coordinate 

Because the cable can stretch and because the amount of 
stretch is a function of the tension T, the independent 
variable s in the above equations, which measures length 
along the cable, is difficult to identify. To reach some 
specific value of s on the cable the observer must know the 
exact state of tension in the cable. A handier variable 
would be some distance between points fixed on the cable 
which was independent of the state of tension. Let us call 
this distance a and measure it on the unstretched cable. This 
new variable is a Lagrangian coordinate in the sense that it 
is fixed to the particles of the cable. For extensible cables 
a simple relation exists between the Eulerian variable s and 
the Lagrangian variable a; that is the stress-strain relation. 
If the stress-strain relation is linear then the relation 
between s and a may be written as 


9s 

9a 



1 + 



/V 



(2.14) 
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where A is the cross-section area of the cable and E the 
modulus of elasticity. The Lagrangian coordinate a may be 
introduced into the dynamic and kinematic equations by a 
change of variables from s to a. There will also be a change 
in cable mass and volume per unit length and a change in 
cable dimensions as well as an effective change in fairing 
mass per unit length and fairing dimensions between the un¬ 
stretched and stretched conditions of the cable. The rela¬ 
tions are 




(2.15) 


where the tilde indicates values for the unstretched cable. 
Making the variable change indicated by equation (2.14) and 
substituting in the relations of equations(2.15) gives for 
the steady configuration equations (2.10a,b) 



for the dynamic time-varying equations (2.11a,b) and (2.9c) 
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:| 


-fh(|^ -v*||-) - (m-pV)g sin(t)*$-T|^*-T*||- -ypCj^c 
[ { (a^+a 2 sin<J>*) 2sin^<t>*+ (aj^+ 2 a 2 sin(})*) cos^(()*}U^u 
+a^sinil)*cos<}t*U^vJ (1+^)^-ipir6^|^ -^pC^c(aj^+a2sin4-*) 

sin4,.U^ (l+2i)-i jfe T ' 0 <2-17^) 

/s 

-fn(|^ +u*|^) - (m-pV) g cos4)*$+|^ -^pCj^c [bj^sin<()*cos(})*U^fl 
+ {2 (b,+b^cos(t)*) cos^4)*+(b-+2b-cos<})*) sin^<|)* }U v] (1+—)^ 

1 Z 1 Z O ATj 

-ipTrB^u*||—^pCjjC(bj^+b2COS(j)*)cos4)*U^(l+^)“^ 2 Xe ^ ® 

(2.17b) 

-fn(|^ +v*|^)+(fh-pV)g<psin<t>*+T*||- -^p7rc^|^ 

-:^pCj^c(aj^+a 2 sin<|>*)U^(l+^)^ w = 0 (2.17c) 

The last term in each of equations (2.17a) and (2.17b) came 
from the time-varying part of the tension in the substitution 
of c for c in the steady equations. Continuing the variable 
change the kinematic equations (2.13a,b) and (2.12c) give 


+ u*3i + (i ii* = 1-11. 
do do do AE at 


M ^ ii* 

aa aa 


_v*|l = - ii 


aa 


at 


a^ 

aa 




T* 


rp* 


(2.18a) 

(2.18b) 

(2.18c) 
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completing the transformation of the equations of motion 
from the Eulerian to the Lagrangian cable coordinate. 

Non-dimensionalization 

The following definitions are made to facilitate non- 
dimensionalization of the cable equations. 

$ = (})* w = w/U^ C = cr/5< n = TPirc^/m 

<() = $ sin $ = u*/U^ T = tg/U^ D = jp Cj^cU^/mg 


6 = U^/lg K = Xg/c 


tp = ij) cos $ = 

u = d/U^ N = T*/mg«, 


d = pV/fn h = m/m 

y = ^pirb^/fh y = AE/fhgil 


V = v/U^ n = T/mgZ 


Substituting these definitions into the equations of the 
previous section gives for the steady configuration 


(l-d)cos'J> - D (aj^+a 2 sin$) sin^Clt—) 




(2.19a) 



(2.19b) 


for the time-varying dynamic equations 



[^(aj^+a2sin4>)sin$ (1+^) ^ 



(2.20a) 
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^ + (1+y) sin^-r^ - 1^ + D(l+^)^ [b-sin$cos$u + 
oT oTo^ Y -L 

{2(b^+b 2 COs$) cos^$ + (b^+ 2 b 2 COs$) sin^$}v] + (1 

+ ^ (b,+b„cos$)cos$ (1+—) ^ n = 0 
2y 1 2 Y 

(l+ri)-|^ + (l+4hC) cos$|^ - n|^ - (1-d) sin$ij^+ 

D(a,+a_sin$) (1+—)^ w = 0 
X 2 Y 

and finally for the kinematic equations 


9? 


cost H + i 


(1+S, 


act) dij) 
3t ~ dc 


V 


0 


lx 

9? 


+ 


sin$ 


a4) 1 an d$ 

9^ Y^ 9t d^ 


0 


aw 

9? 


+ 


cos$ 


liii _ lit _ 

9? 6 Y 9t 


. d$ 
sin4> ^ c|) 


0 


In order to ease the algebraic load the following non 
dimensional drag coefficients are defined and will be 
in later manipulations. 


= D[(a^+a 2 sin$) 2sin^$ + (a^+ 2 a 2 sin$)cos^$] 
= Daj^sin$cos$ 

= D (aj^+a 2 sin$) sin'!’ 


-d) cos$()) 
(2.20b) 

(2.20c) 

(2.20d) 

(2.20e) 

(2.20f) 

used 
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D. = Db, sinfcos'I' 

4 1 

Dg = D [2 {b^+b2COs4') 003^4'+(bj^+2b2COs4') sin^^] 

Dg = D (b^+b2Cos4>) cos4> 

= D (a^+a2sin4') 

Equations (2.19a,b) can be integrated to provide the 
information from which the coefficients of equations (2.20a-f) 
may be determined. Equations (2.20a-f) may be solved by 
techniques described later. Before we can proceed'to a 
solution of these equations, conditions for the behavior of 
the system at the boundaries must be described. 

Boundary Conditions 

The six first order PDE's which comprise the model for 
the cable dynamics require six boundary conditions for their 
solution. Analysis in the characteristic plane shows that 
three conditions at each end of the cable need to be specified. 
Boundary conditions can be either dynamic or kinematic meaning 
that either the forces acting on or the motions of an end of 
the cable may be specified. Different combinations suit 
different problems. For example, a common towing situation 
involves a ship towing a cable and towed body which are small 
compared to the ship. Here there is little feedback from the 
cable dynamics to the ship motions so kinematic B.C.'s are 
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appropriate at the upper boundary of the cable. The motions 
of the towed body are not predictable in advance but the 
forces on the cables resulting from the motions of the towed 
body are, so dyncimic conditions are appropriate at the lower 
boundary. Such a set of B.C.'s will be used for the remainder 
of this study. However, there are other interesting cases, 
notably the opposite of the above case which happens when a 
submerged submarine tows a surfaced buoy in which case the 
kinematic conditions are applied to the lower end of the 
cable and the dynamic conditions to the upper end. More 
complex B.C.'s will result for situations in which the towed 
body develops forces on the tow ship comparable to those 
resulting from its being in a seaway. 

Specifying three dynamic B.C.'s amounts to specifying 
the three components of the force vector applied to the cable 
from the body in each of the three directions of the cable 
coordinate system. Similarly three kinematic conditions 
correspond to motions in each of the three directions. The 
simplest hydrodynamic shape which could be towed would be a 
sphere towed from its center of gravity. Such a shape would 
result in three simple equations which are developed in 
Appendix II. These equations for the dynamic boundary 
conditions at the towed body end of the cable are 

)-^pC^AgU2(^ + ^ ^ (u*^u+u*v*v)) 

o o o 

+ (Mg-pVg)g cos())* -(Mg-pVg)g sin(j)*$ + f^ = 0 (2.21a) 


“^*lt 
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- (M 




o o o 


-(M_-pV_)g sin<l)*-(M_-pV_)g cos())*$+T*+T+f = 0 

l5 15 15 15 y 


(2.21b) 



where M is the mass of the body, its volume, its 
15 15 15 

added mass, C_ its drag coefficient, A„ the cross sectional 

1 15 



to the body from the surroundings which will vanish for most 
circumstances. 

The simplest conditions under which a body might be 
towed would be straight ahead on a calm sea. A slightly 
more complicated and more interesting condition is for the 
tow point to be moving ahead with velocity but with super¬ 
posed sinusoidal oscillations in velocity resulting from ship 
motions. Huffman (1969) considered a towed sphere and straight 
ahead motion of the tow point without oscillation. Unsteady 
motion was introduced by the disturbing forces on the sphere. 
Whicker (1957) considered the case for sinusoidal oscillations 
of the tow point. The former set of B.C.'s are particularly 
appropriate when using the method of characteristics to solve 
the cable dyncimic equations. The latter approach is better 
for solution in the frequency domain. Once the frequency 
response of the tow system is known, the effect of a random 
excitation from the sea surface on the motions of the towed 
body may be evaluated. 
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Suppose the tow point is moving with the following 
velocities in the x^, y^, coordinate system 

u = + 4 a 

o o 

V = h'k 
o 

W = 4 WJ 

o 

where 4lt, and 4t^ are surge, heave, and sway perturbation 

velocities. These perturbation velocities are determined by 
resolution of the six degrees of freedom motions of the ship 
into translational velocities of the tow point. Application 
of the transformation (2.2) gives 

u* + d =. sin(|)*U^ + cos(})*U^$ + sin<|)*4'a - cos<})*/r( 2 . (2.22a) 

V* + d = cos<))*U^ - sin(()*UQ$ + cos(()*4li + sin4>*h'e. (2.22b) 

v) = - i|)cos())*U^ + 4m; (2.22c) 

If the surface perturbations vanish, the result is Huffman's 
B.C.'s and if they are sinusoidal with frequency dependent 
amplitude, the result is Whicker's B.C.'s. Of course, the 
perturbations could be arbitrary functions of time if desired. 
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Equations (2.21) and (2.22) are a simple set of 
boundary conditions. More complicated dynamic conditions 
for the towed body are discussed in Appendix II. Equations 
(2.21) and (2.22) now need to be non-dimensionalized. The 
following definitions are made 

M = Mg/ini!,, k = yg/mil, C = ■^pC^AgU^/mgi!,, a = pVg/fn5, 

X = V = f^/mgJl, Z = f^/mqZ 

Applying the above definitions and separating the equations 
into steady and unsteady terms gives the steady B.C.'s at 
the towed body, that is for C = 0 

$(0) = tan (^) (2.23a) 

N(0) = C /l+(?^)2 (2.23b) 

and the time varying B.C.'s for C = 0 

- (M+k) (^ -cos$^) “C (u+sin^ $u+sin$cos$v) - (M-a) sin4>d>+X = 0 
dT dx 

(2.24a) 

“(M+k) (^ +sin$i^)-C (v+cos^'J’v+sin^cos'J’u) - (M-a) cos'J><j)+n+)'=0 
uT dx 

(2.24b) 

-(M+k) (^ +cos4>^) -Cw+(M-a) sin$4) + Z = 0 (2.24c) 
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The steady kinematic B.C.'s at ? = 1 are trivially satisfied. 
The time varying conditions are 

u = cos^cf) + a (2.25a) 

V = -sin$4) + 3 (2.25b) 

w = -cos<i>i|j + e (2.25c) 

with a = (sin(j)*4'a - cos4>*li'z) /U^, 3 = (cos(})*.4'u + sin(J)*(i'e)/U^, 

e = 4W/U^. 

The problem is now formulated in sets of equations for 
the steady configuration, for the cable dynamics and kinematics 
and for the dynamic and kinematic boundary conditions. 
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Chapter 3 
SOLUTION 


Frequency Domain Analysis 

Suppose the towed system is excited by sinusoidal 
oscillations of the tow point so that in equations (2.25a-c) 
we have for a, S, and e 


a = Re a e^'^^ 

(3.1a) 

3 = Re 6 

(3.1b) 

e = Re F 

(3.1c) 

where Re indicates that the real part is to be taken. 

This 


notation will hereafter be suppressed and it will be under¬ 
stood that the operation of taking the real part must be 


performed in interpreting results. Also t = ^ so v 

o 

wU 

_ _o 

g 


where w is the usual angular frequency. Given the above 
excitation we assume solutions to equations (2.20a-f) of the 
form 

u = u(c)e^'^^ etc. (3.2) 

These solutions are then substituted into equations (2.20a-f) 
resulting in a set of ordinary first order linear differential 
equations in the variables u(c) etc. These equations cast 
into normal form are 
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cos^ (- (l+)j) ivu+cos4>iv4)-D^ (1+^) ^u-D 2 (1+^) ^v-(1-d) 
sinH-[^(l+y)"^+||-]n)- i(l+|)iv^+ g| v (3.3a) 

_ sin«I> (_ (i+p) ivu+cos$iv^-Dj^ (1+^) ^u-D 2 (1+^) ^v-(1-d) 
sin$^-[j|-(l+^)”^+||-]n)+ ^ ivn - H u (3.3b) 


^ = ^(-(l+y)ivu+cos'i>iv(|)-Dj^ (1 +^)^u-D 2 (l+^)^v-(l-d) 
sin$^-[^(l+^)"^+ H-ln) (3.3c) 


^ = ivv + (1+y) sin$iv<j)+D^ (1+^) ^u+Dg (1+^) +(l-d) 

cos $4) + 7 ^ (1+—) ^ n 
2y y 


(3.3d) 


dw _ -cos$ 
dC N 

1 N — 

+ f(l+|)iv4- 


((1+n) ivw + (l+4h5) cosOiviJj- (1-d) sin$iJj+D^ 

+ sinO^ rif (3.3e) 


^ = ^((l+ri)ivw + (l+4h5) cos$ivi|; - (1-d) sin$(}) +D^(l+y)^ w) 

(3.3f) 

Equations (3.3a-f) consist of real and imaginary parts. They 
may be further simplified by dividing the solutions into real 
and imaginary parts. That is 


u 


Uj^ + iUj etc. 


42 




































































The above relations may be substituted into equations (3.3a-f) 
giving a set of twelve equations in the real and imaginary 
parts of the original six dependent variables. The boundary 
conditions may be similarly treated. The set of twelve 
equations with associated boundary conditions divide into two 
sets of coupled equations, one set of eight for the longitud¬ 
inal-transverse vibrations and one set of four for the lateral 
vibrations. These two sets may be written in matrix form with 
their associated B.C.'s as shown below. 
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dc 


dvi 
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‘^4 >r 




d(|)j 


d? 




d^ 


drij 


d^ 



£221 D, (!+-)*■ 

N 1 Y 


cos$ 
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(l-Hi )v 


f COsc^ ^ n+^)^- ^1 


cos^ 
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(1+M )V 


N,r 


. £ost„i,i^a) 


tiiai D,(i+S)^ 

^ N l' Y 




sin^ 
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(1+y)v 


sin$ 


N 


{1+y)V 


N Y 


- i (l+li)v 




siBi D2(i*^)i 


^"N“ D3_(1+y) 


de^ 0 


^(1+P)v 




D, M 1 
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N ^ Y 
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cos4> {1-d) sin$ 

N 




sinll (i-d) 

N ^ ' 


sinj>cos^v 

N 


(1-d) sin^ 

N 


cos j>v 

N 


(1-d)cos$ 




£2|l(^d+!i)-^ + ^1 

N 2y Y 


cos4> (1-d) sin4> 

N 


sin4>cos^v 

N 


sin^ j> (1-d) 

N 


COS't'V 

N 


(1-d) sin4> 

N 


siair^n+Nx-i- 
N 4y' Y 


y6 


V 


lr^(l+N)-2 + 

n‘2y' y’ '3^ 


-(l+lj)sin$v D 


6 

2y Y 


^■27'^^y->'^ ^ t’ 


y6 


" i?' 


it!3,i^.N,-i 4, 

n'2y y’ '3c 


2f'-“7>'' 




u. 


V, 


(3-4) 






(l+y)sin$v 


{l-d)cos$ 


0 




























































-C (l+sin^$) 
- (M+k)v 
-Csin$cos$ 

0 

0 

0 

0 

I 0 

0 0 

0 0 

0 0 

0 0 

1 0 

0 1 

0 0 

0 0 


(M+k)v 
-C(l+sin^$) 
0 

-Csin$cos$ 

0 

0 

0 

0 


-Csin$cos$ 

0 

-C(l+cos^$) 
-(M+k)v 

0 

0 

0 

0 

0 0 ^ 
0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


I 

0 

-Csin$cos$ 
(M+k) V 
-C(l+cos^$) 
0 

0 

0 

0 

Vr(1) 

Vj(l) 

^r(I) 

♦ l(l) 
n^d) 
"l(l> 


-(M-a)sin$ 
(M+k)cos$v 
-(M-a)cos$ 

-(M+k)sin$v 

0 

0 

0 

0 

' 0 

0 

0 

0 

ap^(v) 

aj(v) 

8r(v) 

i 


0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 -cos$ 0 

000 -cos$ 
1 0 +sin$ 0 

0 10 +sin4) 


)- 
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-(M+k)cos$v 

0 

0 ^ 

f 

Ur(0)' 

-(M-a)sin$ 

0 

0 


^j(O) 

+(M+k)sin$v 

1 

0 


7r(0) 

-(M-a)cos$ 

0 

1 


Vj(0) 

0 

0 

0 


*r(0) 

0 

0 

0 



0 

0 
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0 

0 

0 

> 
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Equations (3.4) and (3.5) and equations (3.6) and (3.7) 
for the spatial dependence of motion variables are cast into 
a normal form for linear boundary value problems. 

The terms of the coefficient matrices of these cable 
dynamic equations have a complicated spatial dependence 
arising from the solution of the non-linear steady configur¬ 
ation equations. Therefore a numerical integration scheme is 
essential for their solution. A two part scheme is used. 
First, the steady configuration equations are integrated. 

Then, those results are used in the integration of the 
dynamic equations. The numerical scheme was implemented on 
the IBM 370/155 at the M.I.T. Information Processing Center. 

Steady Configuration Solution 

The equations for the steady configuration (2.19a,b) 
subject to the initial values at ? = 0 given by equations 
(2.23a,b) may be integrated by any of the standard integration 
schemes for a system of non-linear, first order, ordinary 
differential equations. Hamming's Predictor-Corrector method 
was used as provided by IBM (1970) in their Scientific Sub¬ 
routine Package. The problem is straightforward and the 
choice of a particular method of integration is not crucial. 
The results of the steady configuration must be put into the 
coefficient matrix of the dynamic equations when solving them. 
The simplest technique is to make up a table of the values of 
the solution variables, tension and trail angle, as they are 
generated at each step. The table may be accessed with a 
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value of the independent variable. Later when values are 
needed to compute the terms in the coefficient matrices of 
the dynamic equations the table may be entered with a value 
of c and linear interpolation between adjacent tabular values 
of ^ and N used to generate the correct values. A step size 
should be selected which will provide accuracy consistent with 
the requirements of later steps in the computational scheme. 

As mentioned before the results of the integration of 
the steady configuration equations have been considered in 
detail by other investigators, nevertheless, some observations 
are appropriate to assist in understanding the results of the 
investigation of the dynamic equations. Perhaps the most 
striking feature of the steady configuration of a towed cable 
is the sweeping curve it makes in space. Near the towed body, 
curvature is likely to be pronounced with the curve flattening 
as it moves nearer the tow point. Imagine a cable with no 
body at the end. Such a cable will have no curvature. The 
angle at which it tows, called the critical angle, will be 
found by setting d$/d5 = 0 in equation (2.19a) giving the 
transcendental equation 

(l-d)cos$ = D(a,+a«sinO ) sin'!' (1+^)^ (3.8) 

C X ^ C C Y 
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Equation (3.8) shows that the critical angle is the angle 
achieved by the cable as a result of balancing its weight 
with viscous drag forces. Equation (2.19b) may still be 
integrated to give tension along the cable. The initial 
value of tension at ? = 0 will be zero in the absence of a 
body. The tension will vary linearly along the cable as long 
as y>>N which is the usual case. 

Now imagine a body moving at the towing velocity. The 
resultant force from the body's weight and drag will also 
make an angle with the horizontal in general different from 
and certainly independent of the critical angle described 
above. This angle is found in equation (2.23a). When the 
body is attached to the cable the resultant force on the 
body from weight and drag must be matched by the cable which 
is then deformed from its linear trajectory. The effect of 
this body constraint is propagated up the cable toward the 
towpoint dying out as it goes. Therefore the greatest curva¬ 
tures are found closest to the body. Equation (3.8), not to 
mention common sense, shows that the greater the drag the 
smaller the critical angle and therefore the shallower the 
tow for a given length of cable. Similar reasoning applies 
to the towed body. 

In most towed body situations designers seek to minimize 
the length of cable necessary to attain a desired towing 
depth for obvious economic reasons. Therefore a premium is 
placed on drag reduction both of the cable fairing and the 
towed body. This is the principal reason for moving to the 
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complication of a fairing on the tow cable. The other 
important reason is to minimize variable drag and side force 
by minimizing vortex shedding effects. An additional tech¬ 
nique frequently employed to drive a towed body deeper is to 
equip it with a depressor which is a hydrofoil designed to 
develop large downward directed lift forces. A typical steady 
configuration is shown in figure 3-1. When d$/d? ^ 0 due to 
the presence of a body at the end of the cable, the left hand 
side of equation (2.19b) is no longer constant and therefore 
N no longer varies linearly up the cable. 

Dynamics Solution 

Equations (2.29) and (2.31) subject to the boundary 
values expressed in equations (2.30) and (2.32) may be 
integrated similarly to the equations for the steady configura¬ 
tions by any of several standard integration schemes. The 
process is a bit more complicated because boundary conditions 
must be met at both ends of the integration interval instead 
of just at the beginning. Kerney (1971) solves his equations 
by a shooting technique, in which initial values are specified 
as the boundary values at one end, say ? = 0, with missing 
values supplied arbitiarily. Then the integration is carried 
out and the computed values at ? = 1 are compared with the 
desired values given by the boundary conditions. The process 
is repeated with different arbitrary values. Then advantage 
is taken of the linearity of the problem to compute the correct 
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Cable Parameters 


d=0.175 

Y=20.8 

M=4 

D=9.48 

6=0.0066 

b=0 

k=1.035 
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• 
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14 C=0.56 a=2.07 

555 c=2.22 n=0.214 

5 h=0.5 



Figure 3.1 

Typical Steady Configuration 
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missing values. A more direct approach is provided in the 
IBM Scientific Subroutine LBVP (an acronym for Linear Bound¬ 
ary Value Problem). This routine uses solutions of the ad¬ 
joint problem to systematically transfer the boundary values 
at one boundary along the cable to the other end and then 
integrates the resulting initial value problem by a standard 
routine, again the Hamming Predictor-Corrector algorithm, 
HPCG. 


The lateral vibration, cable dynamics equation (3.6) with 
associated B.C.'s in equation (3.7) may be used to show how 
the adjoint technique is applied in LBVP. Equation (3.6) is 
written as 


dw 

dC 


= A w 


where A is the coefficient matrix and w = 


'^R 


wT 



W2 

'f'R 


W3 

KJ 


W4 

. t 


(3.9) 


The adjoint problems are 


dW .Trr 

= -A w 


(3.10) 


with W = 


Wi 

W2 

W 3 

W4 


subject to four sets of B.C.'s 


= Ej^, k=l,2,3,4 such that for k=l, 


^1' 


and for each 
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successive value of k a corresponding adjoint variable boundary 
value at c = 1 is unity, all others vanishing. There are now 
five problems being considered: the original boundary value 
problem and four adjoint initial value problems, one for each 
value of k. The adjoint problems may now be integrated by 
HPCG from ? = 1 to ? = 0, thus obtaining values for the ad¬ 
joint variables at C = 0 in each of the adjoint problems, 
d Ic d 

Now w^] = w] where the superscript k refers 

th 

to the k^ adjoint problem. Carrying out the differentiation 

on the R.H.S. and substituting from equations (3.9) and (3.10) 
d k 

gives Wj^] = 0 which may be integrated giving 


EW^w.l , - EW^w.l 

1 X 'C=1 . 1 1'C=0 


ecus. The k^ adjoint B.C 


^ ^i -ilr=0 = ^klr=l 


0 since equation (3.9) is homogen- 
at c; = 1 may be applied giving 

(3.11) 


There four of these equations, one for each value of k. Now 
the B.C.'s from the original problem are used. Suppose they 
are written as 


Bw|^^q = 0 (3.12a) 

Cwl^^l = F (3.12b) 

where B and C are 2x4 coefficient matrices, that is (3.12a) 
and (3.12b) give two equations each. Substituting from 
(3.11) into (3.12b) for each value of k gives two equations 
now in terms of the unknowns w|^_q which combined with 
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equation (3.12a) gives four simultaneous equations in four 
unknowns, namely the initial values of the original problem. 
The coefficients of the terms in these equations are all 
known being combinations of the towed body parameters, the 
trail angle at c = 0 and 1 obtained from the steady con¬ 
figuration solution, the frequency of interest and finally 
the solutions of the adjoint problems. These four equations 
may now be solved simultaneously to give the initial values 
for the integration of the original equations by HPCG. The 
transverse-longitudinal problem is handled similarly. 

Solution Difficulties 

The numerical scheme presented above is useful over a 

wide range of practical interest. It suffers, however, from 

a purely computational malady as frequency, cable length or 

viscous damping is increased. The character of the solution 

of equations (3.4) and (3.6) is determined by the eigenvalues 

of the coefficient matrices. Were the coefficients constant 

instead of space varying, the solution to the equations would 

be simply found by calculating the eigenvalues and their 

associated eigenvectors. Solutions would be of the form 

u(c) = E c.e^^^ where c. are the eigenvectors and the 
i 

eigenvalues. Such a procedure is in fact used in Appendix III 
to verify the validity of the numerical method in certain 
special cases wherein the coefficient matrices of equations 
(3.4) and (3.6) are constant. The coefficient matrices being 
space varying means that their eigenvalues change along the 
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cable but that does not affect the essence of this argument. 
Suppose an eigenvalue, has a large positive real part. 

Then the quantity, e^^, occurring when the B.C.'s at ^ = 1 
are satisfied, has a very large positive real part. Therefore 
arbitrary constant in c^ must be made extremely small in order 
to match a reasonable boundary value. 

Physically, the solution consists of dying and growing 
waves emanating from the excitation at 5=1. Since the towed 
body end of the cable is free to move about, subject to the 
constraint of the dynamic B.C. at 5 = 0, a little of the 
growing solution is needed to satisfy this B.C. The growing 
solution eigenvalue has a positive real part and the dying 
solution a negative real part. As the solution is more and 
more damped, because the frequency, steady tow velocity, or 
damping coefficients have increased, or because the distributed 
damping along the cable has had a greater effect because the 
cable is longer, less and less of the growing solution is 
needed to satisfy the B. C. at5=0. If so little of the 
solution is needed that the value of the arbitrary constant 
in the eigenvector cannot be made small enough because of 
round-off in the computer, the numerical scheme will fail to 
provide a solution. Kerney (1971) apparently had similar 
difficulty. Practically, this difficulty seems not to be 
egregious since the more highly damped situations are exactly 
those in which the least motion is transmitted to the towed 
body. 
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In the computional scheme employed using the subroutine 
LBVP the difficulty manifests itself in repeated attempts to 
subdivide the integration step size in order to compute the 
solution within a small error bound. The error bound is small 
in order to test appropriately for loss of significance in the 
Gaussian elimination algorithm used to solve the set of 
simultaneous algebraic equations generated with the adjoint 
solutions for the boundary values at ? = 0. 

Reducing round-off error by moving to higher precision 
is a way of enlarging the range of parameters for which solu¬ 
tions may be found using the method described in this thesis. 
Comparison of solution efforts using single and double pre¬ 
cision shows that solutions which fail in single precision go 
forward in double precision. FORTRAN compilers are becoming 
available in extended precision so the range may be increased 
still further without resorting to the complexity of multiple 
precision programming. Table 3.1 shows the value of the 
largest positive real part of the eigenvalues of the coeffici¬ 
ent matrix computed for the point on the cable where the 
largest values occur when the solution method fails in single 
and double precision. Table 3.1 shows that the solution 
proceeds for much larger positive real parts of the eigen¬ 
values corresponding to a larger range of problem parameters, 
in this case frequency, when it is executed in double 
precision. 
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CABLE d=0.125 y= 2775.0 M=2.755 C=0.373 a=1.38 D=9.48 6=0.00438 

DATA b=0.555 c=2.22 y=0.214 k=0.69 ri=3.44 ^=0.5 h=0.5 

precision 

step size 

error requirement 

positive real part 
of eigenvalue at 
solution failure 

single 

0.005 

0.00001 

13.43 at v=0.175 

double 

0.004 

0.000004 

20.90 at v=0.375 


Comparison of Single and Double Precision Programs 

Table 3.1 

While the solution difficulty described above is 
troublesome, the method of solution is useful since many 
practical situations may still be considered. Furthermore, 
the difficulty can be put off to some degree by increasingly 
precise computation. 

Future efforts directed toward more sophisticated solution 
schemes which avoid the above difficulty are merited. When the 
excitation frequency is greatly increased, the short wave 
length approximation will apply. That is the steady configura¬ 
tion trail angle and tension will change slowly over a wave 
length so the approximation is made that they are constant 
over a range of several wave lengths. The WKBJ method 
(cf. Morse and Feshbach (1953)) then applies. Using the same 
approximation Milgram (1971) suggests that a Green's function 
approach will work. Another useful line of investigation 
might be to look for ways to control the growing exponential 
solution in the numerical scheme, perhaps by segmenting the 
integration range in some fashion. 
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Chapter 4 
RESULTS 

The utility of the formulation and solution scheme for 
towing cable dynamics presented above can be demonstrated by 
applying it to a plausible cable design and examining the 
results. First a standard reference spectrum of tow point 
motions is prescribed. For the results described in this 
chapter the towpoint translational perturbation velocities in 
the y^, coordinate system are specified. The 

amplitudes of these perturbation velocities are unity when 
normalized with respect to the tow velocity U^. A typical 
input perturbation is Uq/Uq = = 1 cosvx. Because the 

system described is linear the results for any arbitrary 
spectrum may be readily obtained from the results for this 
flat spectrum with a normalized amplitude of unity. Second, 
this disturbance is applied to a basic towing cable and body 
design which is representative of real towing situations. 
Appendix IV contains a description of this system and shows 
how the non-dimensional parameters were obtained. The basic 
trial cable design envisions towing a 40 in. diameter sphere 
at 16 ft. per sec. with a 1200 ft. steel cable. 

The greatest interest in these matters is in the motion 
of the towed body and in the largest unsteady values of 
tension in the cable. Such values are calculated and presented 
graphically. However, motions and tension anywhere along the 
cable could be determined as well. The amplitude of the 
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motions is of more consequence than phase so only amplitudes 
are presented, again normalized with respect to steady tow 
velocity for the resulting velocities or with respect to 
the weight in air of the entire cable for the tensions. 
Furthermore, the output velocities have again been trans¬ 
formed from the cable coordinate system x, y, z back to the 
locally vertical system x^, y^, z^ via the transformation in 
equation (2.2). 

After the motions and tensions of the basic cable are 
established, significant parameters are varied over a wide 
range. Cable lengths from 300 to 1800 ft., steady tow 
velocities from 0.1 to 48 ft. per sec., and moduli of 
elasticity from 2.5 x 10** psi to virtually inextensible are 
considered. Variations were made one parameter at a time 
using the basic cable as a reference. Other parameters 
could also be varied and indeed probably would be in a 
design situation. Length and velocity variations seem most 
significant since for any design these parameters are likely 
to vary during the use of the system. Modulus of elasticity 
was varied because the predicted results of moving to 
springier cables are particularly interesting. Finally, 
researchers have frequently ignored the added mass of the 
cable in previous investigations so the basic cable was 
considered sans cable added mass and results compared. 

A physical feel for a complicated problem like this may 
be elusive. But the idea of a damped oscillator is basic to 
the physics of this interaction and so the natural frequencies 
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are of interest. Bearing in mind that this idea is not 
being applied as a tool of analysis but rather as a way of 
understanding, natural frequencies of various idealisations 
of the problem are calculated, compared to sample results 
and discussed. 

Basic Cable System 

The results from applying the solution method to the 
basic system are shown in figure 4.1. The magnitude of the 
unsteady tension n^, at the tow point where it is greatest 
rises approximately quadratically with increasing frequency 
in the higher frequency range. This is expected for the 
nearly inextensible steel cable and reflects the increasing 
inertial forces at higher frequencies. The magnitude is 
largest at the tow point because the cable there must carry 
the load of its entire length and the body as well. The 
amplitude of lateral velocity w^^ decays rapidly with in¬ 
creasing frequency showing the effect of heavy viscous damping, 
v^j^ and u^j^ cannot decay substantially because they are in 
the plane of the steady configuration and the near inexten¬ 
sibility of the cable implies that most of the tow-point 
motion will be transferred to the body. Some energy is re¬ 
moved from the cable by the viscous damping but most is 
transmitted down the cable to the body. The relation between 
the amplitude of the vertical and horizontal components seems 
to be a strong function of the relation between the trail 
angle at the tow-point and towed body, as will be apparent 
in later figures. 
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Cable Parameters 


d=0.125 y=4160.0 

D=9.48 6=0.0066 

k=1.0315 ri=3.44 

normalized velocity 
magnitude 


M=4.14 C=0.56 
b=0.555 c=2.22 
5=0.5 . h=0.5 


a=2.07 

y=0.214 



Motions and Tensions - Basic Cable 
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Variation of Cable Parameters 


The results of varying the length of the cable are shown 
in figures 4.2 and 4.3. Table 4.1 shows the values of the 
cable parameters used. As cable length increases, motions 
are diminished because there is more total damping, that is 
the same drag per unit length is acting over a longer length 
of cable. The vertical and horizontal motions tend toward a 
limit defined by the relation of the angle of the resultant 
force at the towed body and the cable critical angle. As the 
cable becomes longer and longer the angle at the towpoint 
approaches the critical angle, and so the motions approach a 
limiting case. As length increases the lateral motion of the 
towed body becomes more and more rapidly damped with increasing 
frequency. 

The variations of steady tow velocity give more dramatic 
results. Table 4.2 and figures 4.4, 4.5, and 4.6 show the 
results of varying velocity from 0.1 ft. per sec. where the 
cable is nearly hanging vertically below the tow point to 
48 ft. per sec. where it is towing nearly horizontally. 

The two striking features are the appearance of resonances 
when the linearized viscous damping is quite small correspond¬ 
ing to being small and the way in which the predominant 
motion shifts from vertical when the cable is hanging vertically 
to horizontal when the cable is extended horizontally. In 
both cases the vibrations normal to the cable are being 
damped while the longitudinal vibrations are transferred 

directly to the towed body. The lateral motion at the 
smallest towing velocity of 0.1 fps may be in error at the 
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1 ft. 

Y M 

C a 

6 


k 

300 

16640.0 16.56 

2.240 8.28 

0.0264 

4 

.14 

600 

8320.0 8.28 

1.120 4.14 

0.0132 

2 

.07 

900 

5550.0 5.52 

0.746 2.76 

0.00879 

1 

.38 

1200 

4160.0 4.14 

0.560 2.07 

0.0066 

1 

.03 

1500 

3330.0 3.31 

0.448 1.656 

0.00528 

0 

.828 

1800 

2775.0 2.755 

0.373 1.380 

0.00438 

0 

.690 

Changes from Basic Cable 

Parameters for 

Length Variations 



Table 4.1 




U f .p 
o ^ 

. s. C 

D 

6 



0.1 

0.0000219 

0.00037 

0.00000026 


1 

0.00219 

0.037 

0.000026 



3 

0.01967 

0.333 

0.000235 



6 

0.0786 

1.331 

0.000938 



12 

0.314 

5.33 

0.00375 



18 

0.706 

11.92 

0.00843 



24 

1.256 

21.30 

0.0150 



48 

5.020 

85.30 

0.0599 




changes from Basic Cable Parameters for Tow Velocity Variations 

Table 4.2 


E psi X 10® 

Y 

CO 

99999999.0 

■ 10 

2080.0 

5 

1040.0 

0.5 

104.0 

0.1 

20.8 

0.05 

10.4 

0.025 

5.2 


Changes from Basic Cable Parameters for Elasticity Variations 

Table 4.3 
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Lateral Motion and Tension - Length Variation 
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Cable Parameters - see table 4.2 


normalized 
velocity magnitude 



Figure 4.4 

Horizontal Motions - Steady Velocity Variation 
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Cable Parameters - see table 4.2 

normalized 



Figure 4.5 

Vertical Motions - Steady Velocity Variation 
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normalized Cable Parameters - see table 4.2 



Figure 4.6 

Lateral Motions - Steady Velocity Variation 
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higher frequencies because the low reduced frequency criterion 
in the derivation of the lateral hydrodynamic forces is being 
violated. Note that the shapes of curves near resonances in 
the various figures are only approximate due to the narrow 
bandwidths sometimes encountered. 

Most investigators of tow cable dynamics have assumed 
that the tow cable was inextensible. This certainly seems a 
reasonable assumption for a steel cable. Huffman (1970) , how¬ 
ever, recognized that if a time domain solution by the method 
of characteristics is to be attempted the cable must be 
assumed extensible so that the initial value problem is well- 
posed, an unusual case where complicating the problem makes it 
more tractable. The results shown in figure 4.7 are another 
reason why cable extensibility is important. The more elastic 
the cable the less it couples the surface motions to body. 
Potential resonances are unlikely to be a significant problem 
because there is considerable viscous damping present. In 
figure 4.7 vertical and horizontal normalized velocities of 
the towed body are shown for cables with moduli ranging from 
2.5 X 10** psi to inextensible and the associated natural 
frequencies for an idealized mass-spring oscillator of 
similar proportions. Associated cable parameters are shown 
in table 4.3. For the more elastic cables there is a great 
deal of coupling between the transverse and longitudinal modes 
of vibration. Transverse damping is much greater than long¬ 
itudinal damping. Energy is coupling from the stress waves 
into the transverse tension waves and then being damped out 
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Cable Parameters - see table 4.3 


normalized 



Figure 4.7 

Horizontal and Vertical Motions - Elasticity Variation 
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because of the viscous effects. Bringing the natural frequen¬ 
cies of the longitudinal stress waves and the transverse tension 
waves into the range of the forcing frequency facilitates damp¬ 
ing and reduces the motions of the towed body. Figure 4.7 
shows how the perturbation velocities of the towed body are 
greatly attenuated for the lower values of modulus of elastic¬ 
ity. This process might be abetted in a real system by inserting 
springs at intervals along the cable or by choosing a cable 
material with a fairly low modulus of elasticity. 

A solution was generated for the basic cable system with¬ 
out added mass. Results were within 1.0% of those of the basic 
system with added mass. Ignoring cable added mass is justified 
for cable systems similar to the basic one studied here. 


Natural Frequencies 

The natural frequencies of the modes of oscillation for 
three idealizations of the cable dynamics problem may be 
readily calculated. The first such idealization is that of a 
compound pendulum. Ignoring viscous effects on the cable and 
considering it to be straight, the linearized pendulum equa¬ 
tion for transverse motions is 

(M+k+-J-(l+y) + (N(O)mgil^ +^mgJl^ sin4>)(j) =0 (4.1) 

giving 

tOp^ = /f and for lateral motions (4.2) 


/N (0)+sin4>/2 /g 

Pil / M+k+i (1+n) / a 
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where the restoring force of the pendulum is composed of 
gravitational forces on cable and body and viscous drag 
forces on the body. A second idealization is to suppose that 
the cable is a straight stretched string with supports at each 
end offering infinite impedance, subject to constant tension 
and again that viscous effects are negligible. Then the 
appropriate equation is the wave equation and the well known 
results (cf. Morse (1948)) for the fundamental natural 
frequencies of transverse and lateral motion are 


0) 


Tt 


TT 


/ N(o) 

/ (1+y) 



and w 


Tfi, 


TT 


/N(o) 

/ (1+n) 



(4.4,5) 


Finally the extensibility of the cable gives rise to 
longitudinal stress waves which cause the cable system to re¬ 
semble a mass spring oscillator, the natural frequencies of 
which may be calculated from the Rayleigh method (cf. 
Timoshenko and Young (1955)) to give 

“s ' /wk+TT /I ^ spring constant K = 

(4.6) 

again ignoring viscous effects on the cable and assuming it 
to be straight. Figures 4.8 and 4.9 show the consequences of 
allowing the fairing drag coefficient of the basic cable to 
vanish. The appropriate natural frequencies are shown for 
comparison. Figure 4.10 shows the horizontal and vertical 
velocities at the towed body for a similar situation except 
with a rather small modulus of elasticity. Comparing figure 
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, . , Basic Cable with = 

normalized D 



Horizontal and Vertical Motion with Cj^ = 0.0 


73 












i * 



I 





0.0 


normalized 


Basic Cable with = 



Lateral Motions and Tension with = 0.0 
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0.0 


Basic Cable with E = 0.05 x 10® psi, = 


normalized 



cOp^ angular frequency w 

Figure 4.10 

Horizontal and Vertical Motions - Very Elastic, = 0.0 
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4.10 with 4.7 shows the effect of heavy viscous damping dis¬ 
tributed along the cable. In the case shown in figure 4.10 
there is viscous damping only on the towed body so modified 
resonances are apparent in the frequency response curves. In 
figure 4.7 for the same modulus of elasticity the damping is 
very large and the resanances have vanished. Even in figure 
4.10 though, the advantaae of reduced cable modulus in 
limiting the motions of the towed body is apparent. 

Figures 4.8, 4.9 and 4.10 indicate that the basic motions 
of the system are those of various simple harmonic oscillators 
coupled in complicated ways. When the fairing drag coefficient 
is restored to a reasonable value the system becomes heavily 
damped and resonant behavior at the natural frequencies 
vanishes. 
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Chapter 5 
CONCLUSION 

Design Impact 

The results of this thesis may be summarized by mention¬ 
ing their possible impact on the design of a towed cable 
system. The first result is a formulation and a method of 
analysis which may be easily implemented, using the FORTRAN 
Scientific Subroutine Package. The formulation is complete 
considering motions in three dimensions with a faired, exten¬ 
sible cable. The method allows the designer to use, directly, 
existing studies of ship motions in the frequency domain. 

Using the formulation and method the designer can find the 
towed body motions for given tow-point motions. The difficulty 
with the method is its susceptibility to overflow on the com¬ 
puter caused by the round-off error interacting with the grow¬ 
ing exponential solution. This difficulty is mitigated by 
three factors: it is associated with increasing viscous damp¬ 
ing and does not foreshadow any reasonance phenomenon, it may 
be put-off by employing more precise computation, and the 
range of application of the method already includes many 
practical designs. A second result is that fairing added 
masses are not important and may be ignored in typical designs 
without appreciable change in resulting motions. Nevertheless 
the formulation of cable fairing hydrodynamics of Appendix I 
is available. A third result is that designers may use fairly 
elastic cables and take advantage of the heavy viscous damping 
to minimize motions by creating a system in which the fundamen- 
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tal natural frequencies of transverse tension waves and 
longitudinal stress waves are small and are in the range of 
excitation by the tow point motions. 

Suggestions for Future Work 

The present need in cable dynamics is for experiments to 
t^st the various studies which have been performed. Experi- 
m ents ar e needed in two areas, cable system dynamics as 
formulated in the body of this thesis and the unsteady hydro¬ 
dynamics of fairings as formulated in Appendix I. Froidevaux 
(1968) has devised a scheme using an inertial platform to 
measure towed body motions which may serve as a useful start. 
The point is that the theory of this thesis and its antece¬ 
dents have not been systematically evaluated by experiment ex¬ 
cept for the steady drag forces on the cables and the steady 
configurations. 

Theoretical work is needed to compare the method of char¬ 
acteristics solutions of Schram and Reyle (1968) and Huffman 
(1970) with the frequency domain solutions of Kerney (1971) 
and of this thesis. A general formulation of the unsteady 
hydrodynamics of cable fairings in three dimensions is needed 
to replace the compound approach of Appendix I. Finally more 
general dynamic boundary conditions need to be developed to 
handle the case of a small tow ship and a large body where 
motions of the body affect motions of the tow ship and cases 
where the towed body is a more general shape than the sphere 
employed in this thesis. 
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Appendix I 

HYDRODYNAMICS OF CABLE FAIRINGS 


Introduction 

No general formulation of the hydrodynamics of cable 
fairings exists. Previous investigators have largely confined 
their attention to the drag on a fairing resulting from a 
steady flow co-planar with the configuration of the cable. 

The problem of interest in these investigations was to find 
the drag forces normal and tangential to the fairing as a 
function of the trail angle, that is the angle of inclination 
of the fairing to the inflow velocity vector. Most investi¬ 
gators have developed loading functions. These functions are 
the normal and tangential drag forces normalized with respect 
to the drag force on the fairing when the trail angle is 90°. 
Casarella and Parsons (1970) present most of the significant 
results in a table of loading functions for faired cables. 

The usual approach has been to conduct experiments and then 
fit a curve through the data. Because the normal and 
tangential drag components are basically related to the 
squares of the normal and tangential velocity components, the 
usual effort is directed at constructing curves from combina¬ 
tions of sines and cosines of the trail angle. Two investi¬ 
gators, Choo (1970) and Calkins (1970), have attempted 
theoretical formulations of the loading functions although 
their methods still require some experimental input. 
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A more comprehensive formulation of cable fairing hydro¬ 
dynamics than discussed above is necessary in an investigation 
of the dynamics of a faired cable in a towed system undergoing 
motion in three dimensions. An understanding of a three- 
dimensional unsteady problem as opposed to a two-dimensional 
steady problem is required. Added mass as well as viscous 
effects and the effects of motions out of as well as in the 
plane of the steady tow configuration must be considered. In 
lieu of a general formulation which may be quite difficult the 
following piecemeal approach has been used to understand the 
hydrodynamics of a cable fairing segment undergoing unsteady 
motions in three dimensions. The contributions of added mass 
to the forces and moments on the fairing segment are written 
following Imlay (1961). The loading function formulation of 
the in-plane viscous drag forces is assumed to apply to the 
unsteady case. The out of plane added mass and viscous forces 
are calculated from thin wing theory following Robinson and 
Laurmann (1952) after imposing some additional restrictions on 
the fairing. These results are then superposed to give forces 
and moments on the cable fairing segment. Unfortunately there 
are no experimental results available with which to test the 
compound theory described above. 

Before proceeding to a detailed discussion of the com¬ 
pound hydrodynamic theory of cable fairing segments some 
simple assumptions are made about the attachment of cable 
fairing segments to the cable. The fairing is assumed to be 
attached to the cable by a frictionless bearing such that 
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forces are communicated to the cable but moments are not. 

For most fairing designs this is an excellent assumption for 
forces normal to the cable, less so for tangential forces. 

In most cases a fairing is somewhat free to slide along the 
cable and is only restrained from doing so by the proximity 
of adjacent segments. In a few designs segments are attached 
to collars fixed to the cable and the above assumption is 
then also correct for tangential forces. In a few cases not 
considered in the following theory fairings are bonded into 
the structural element of the cable so that moments exerted 
on the fairing about the cable axis are directly transmitted 
to the cable. As shown by Abkowitz (1965, 1967) these designs 
are subject to a kind of flutter instability and can apparently 
also lead to remakrable kiting effects. However, in many 
designs, probably to avoid the expense of attaching collars 
to the cable to support each fairing segment separately, 
adjacent segments are free to bunch together or are prevented 
from bunching by a flexible separator between segments. 

Either way a moment-side force coupling might be envisioned 
in so far as a segment is prevented from aligning with the 
local flow. In fact this complicated kind of attachment 
might be modeled as a pair of coupled cable dynamic problems, 
one for the structural cable and one for the fairing treated 
as a different cable with coupling through the transmission 
of normal forces from fairing to structural cable. Such an 
adventure awaits a future investigator. 
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Figure I-l shows a typical fairing segment, associated 
cable and a coordinate system with origin on the centerline 
of the cable. A moderate simplification of algebra results 

I 

later if the cable radius is assumed small compared to the 
chord length of the fairing so that the y axis of the coordinate 
system lies essentially along the leading edge of the fairing. 
The fairing is assumed to be symmetrical about its centerline 
plane. Real cable fairings are designed to be symmetrical 
but sometimes acquire a permanent set from being wound on a 
winch drum, producing a slight asymmetry. Adjacent segments 
are assumed to be infinitesimally far apart so that there are 
no gaps in the fairing. Actually gaps will exist and some 
complicated flow will go on in the gaps. None of the above 
assumptions about the attachment of the fairing to the cable 
is thought to cause significant disparity between theoretical 
and actual performance of the towed cable systems. 

Added Mass 

The added mass tensor appropriate for the arbitrary trans¬ 
lations and rotations of the differential fairing segment in 
the cable coordinate system with its origin at the segment 
mid-plane, as shown in Figure 1-2, is given below. The seg¬ 
ment length ds is large compared to the chord of the fairing 
but small compared to the radius of curvature of the cable. 

The notation of Imlay (1961) is used in which each term is a 
hydrodynamic force or moment derivative per unit length of 
fairing with respect to the subscripted variable. Terms which 
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Figure I-l 

Typical Fairing Segment 



Rotational and Translational Velocities of the Fairing Segment 
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vanish do so either because of the geometrical symmetry of the 
fairing segment about its centerline plane and about its mid¬ 
plane or because the fairing is continuous in the y direction 
of the cable coordinate system. 


0 

0 

0 

0 

0 


0 0 

0 0 

0 Z* 

w 

0 0 

0 M* 

w 

0 0 


0 0 0 

0 0 0 

0 Z‘ 0 

q 

K* 0 0 

p 

0 M* 0 

q 

0 0 N* 

r 


(I-l) 


where X, Y, and Z are forces in the x, y, and z directions and 
K, M, and N are moments about the x, y, and z axes respectively. 
The forces and moments on a rigid, isolated fairing segment of 
length ds retaining only linear terms in the perturbation 
variables are then given by 


= X* ds u 
u 

Y^^^ = X* ds u* B 
u 

Z= Z* ds ^ + Z• ds q - X* ds u* q 
w q u 

(I-2a-f) 

= Z* ds V* w + Z* ds V* q + K* ds p 
w q P ^ 

= Z* ds (w-u*q) + M* ds q - (Z*-X*)ds u* w 
q ^ q w u 

= -X* ds (u*v*+u*'0+v*u) + N* ds ^ 
u r 
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Some simplifications result from considering how the 
above hydrodynamic forces and moments on a fairing segment 
from added mass interact with the cable. Equations (I-2d) 
and (I-2f) give moments about axes normal to the cable element 
and are therefore bending moments. The basic cable assumption 
is that bending moments are negligible. For bending moments 
to contribute to cable deformation the moment at any point in 
the cable must be finite, that is of 0(1). Equations {I-2d) 
and (I-2f) give bending moments of 0(ds) and therefore may be 
ignored in considering the cable dynamics. 

Equations (I-2a) and (I-2b) will be combined with the 
results of considering viscous drag forces to give the com¬ 
ponents of the total hydrodynamic forces on the cable in the 
X and y directions. may be estimated by assuming that the 

added mass is similar to that of an elliptical section of the 
Scime thickness. This assumption gives = -^pirb^ when b 
is the maximum breadth of the fairing. 

Equations (I-2c) and (I-2e) will be used in conjunction 
with the results from thin wing theory to give the contribu¬ 
tion to lateral forces and moments due to added mass caused 
by the thickness of the fairing. These corrections to the 
thin wing theory will be the terms in equations (I-2c) and 
(I-2e) involving X^. 
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Loading Functions 

Viscous drag forces on cables have customarily been 
given by investigators in the form of loading functions. In 
the two-dimensional steady problem as shown earlier in this 
paper the dependent variables are tension and trail angle. 

The trail angle gives the orientations of the cable relative 
to the steady flow. If the cable is oriented normal to the 
flow corresponding to a trail angle of 90° the drag force 
will be given by 

R = -Cp ip c ds Vlv| 

where is a drag coefficient depending on Reynold's nvimber, 
C the chord length of the fairing, ds the span of the fairing 
segment and V the fairing velocity vector. For the equations 
of motion of the cable element the normal and tangential 
components of the drag force are of interest. As the cable 
is inclined to the flow these components will vary as 
functions of the inclination angle, 6. Normal and tangential 
loading functions may then be defined as 


f 

n 


D n (6) 

normal 

-R 


(I-4a) 


D. 4.- 1 (6) 

_ tangential 

-R 


(I-4b) 


where f^ and f^ are the normal and tangential loading functions, 

Vrmal ^tangential normal and tangential compon¬ 
ents of the viscous drag force and R = |r|. f^^ and f^ have 
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been measured experimentally and expressed in terms of sine 
and cosine series. Casarella and Parsons (1970) give the 
forms of f^ and f^ which have been developed by various 
investigators. The theoretical formulations of Calkins (1970) 
and Choo (1970) are of interest but are no more useful in 
the development of the theory of this paper than the experi¬ 
mental results. In the absence of new definitive experimental 
results for a variety of fairing shapes and Reynold's numbers 
the loading functions developed by Whicker (1957) seem as 
relevant as any. These loading functions are 

f = a, sin 0 + a_ sin^ 9 (I-5a) 

n 1 2. 

f^ = bj^ cos 0 + b 2 cos^ 0 (I-5b) 

b b 

where the coefficients are a, = 1 - —, a» = —, b, = 0.386 - 

1 c 2. c L 

0.303 b- = -0.055 + 0.020 -. 
c 2 c 

The following development of the unsteady viscous drag 
forces on a cable fairing segment using Whicker's loading 
functions could be applied to other loading functions as well. 
In the steady flow case the angle of inclination of the cable 
to the flow, 0, is identical to the trail angle, (j). The 
unsteady case is more complicated. First, the assumption is 
made that the drag coefficient for the steady flow applies to 
the unsteady case. Second, the local inflow velocity will 
depend on unsteady as well as steady terms, similarly the 
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I 

angle of inclination of the fairing to the flow will involve 

more than the trail angle alone. Neglecting for the moment 

motion in the direction at any point along the cable the 

local flow will have components in the x^, coordinate 

system ^o' ^o* magnitude of the inflow will then be 

|v| = /U^ + 2U u + + v2. Then as shown in figure 1-3 the 

o o o o o 

angle of inclination to the flow 0 will be given by 

e 4) - 3 ( 1 - 6 ) 

where 3 = and 4) = 4)* + $. Furthermore the coordinate 

transformation from x^, y^, to x, y will give the local flow 
velocities in the cable coordinate system. They are 

u = u* + 0 = (sin4)* + cos4i*$) ~ cos4>*v^ (I-7a) 

V = V* + V = {cos4)* - sin4)*$) (U^+u^) + sin4)*v^ {I-7b) 

Equations (I-7a) and (I-7b) give 

u^ = sin4)*u + cos4>*v (I-8a) 

v^ = -cos4)*u + sin4)*v + (I-Sb) 

Combining equations (I-8a,b), (1-6), (I-5a,b), (I-4a,b) and 

(1-3) will give expressions for the normal and tangential 
drag forces in terms of the dependent variables of the cable 
dynamics formulation. 
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Inflow velocity 
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Figure 1-3 

Inflow Velocity and Components 



Figure 1-4 

Lateral Velocity and Force Diagram 



Figure 1-5 
Fairing Dimensions 









I 


D 


normal 


= -Cp^pc ds [ {aj^+a2sin4)* }sin<})*U^+{ (aj^+a 2 sin<j)*) 


2 sin^4)*+ (a^+2a2sin4>*) cos^^** }U^u+aj^sin<})*cos4)* 



(I-9a) 


D 


tangential 


= c ds [ {bj^+b 2 COs 4 )* }cos<})*U^+bj^sin(})* 


cos(|)*U^d + {2 (bj^+b 2 Cos<()*) cos^4)* + 


(b^+2b2Cos4)*) sin^ (j)* }U^^] 


(I-9b) 


Now that an expression for the normal force has been developed 
the effect of the z direction motion may be included. As shown 
in figure 1-4 the normal component of the viscous drag on the 
fairing has components in both the x and z directions of the 
cable coordinate system. The z component results from the 
small component of velocity in the z direction, v). The normal 
force is along the direction of the normal velocity, therefore 
it will have a z component proportional to the z direction 
velocity component. Similarly the x direction component of 
the normal force will be proportional to the x direction 
velocity component. Retaining only terms which are steady 
or linear in the perturbation quantities gives 


Z 


= -Cjj ^p c ds{aj^+a 2 sin 4 )*}U^W 


(I-lOa) 


90 








= -Cj^ ip c ds [ {aj^+a 2 sin(|)* }sin(|>*U^ + 

{ {aj^+a 2 sin(l>*) 2sin^(l)* + (aj^+ 2 a 2 sin(J)*) cos^(|)*} 

U^u + a^^sincj)* cos(})* U^v] (I-lOb) 

The component of the viscous force in the y direction is just 
the tangential force of equation (I-9b). 


Thin Wing Theory 

An effect of the flow about the fairing in the x-z plane 
is to develop a side force and moment on the fairing due to 
its behavior as a hydrofoil. The fairing is assumed to be a 
thin symmetrical foil moving approximately along its x axis 
with velocity u* but having also a small amplitude motion in 
the lateral or z direction. The resulting flow is assumed to 
be unaffected by the spanwise flow in the y direction. The 
results presented by Robinson and Laurmann (1952) for the two- 
dimensional motion of a thin wing with constant forward speed 
then apply to the fairing. These results include added mass 
and vortex effects. The equations for the hydrodynamic 
derivatives per unit length are 


Z=Z w+Z*w+Z a+Z q+Z*q 
w w a q q 

M = M + + M a + M q + M*q 

w w a q q 


(I-lla,b) 
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9 ot 

where a is the fairing rotation angle and q = Corrections 

must now be supplied for the thickness of the fairing and for 
the effect of rotations about the x and z axes with velocities 
p and r. Examination of equations (I-2c) and (I-2e) show the 
presence of a term involving X^. This term is interpreted as 
a correction for thickness. Equations {I-2c) and (-2e) show 
no added mass effects proportional to p or r, so the only 
contribution from these rotations could be from vortex effects. 
Such contributions are assumed to be of higher than first 
order in perturbation variables and are therefore not considered. 

When considering the interaction of cable motion with the 
cable fairing the usual assumption (cf. Whicker (1957)) is 
that the fairing instantaneously aligns with the local 
direction of flow. In fact, of course, the fairing must 
respond to the dynamic interaction of the hydrodynamic forces 
from thin wing theory, thickness effects, and viscous drag 
with gravitational, inertial and external forces, the external 
force being the lateral force of contact between the cable 
and the fairing. The dynaimic equations for the side force 
and moment per unit length acting on a fairing segment may 
now be written with the parameters of the lateral motion 
a, and <p as dependent variables. The equilibrium of 
forces gives 

-m [w+pv*-qu*-x^q] + (m-pl/)g (acos<{)*+i|isin<{)*)+Z +Z^w+Z^^+ 

Z^a+Zgq+Z^4-X^u*q+Z^ = 0 (1-12) 
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where m is the mass per unit length of fairing, \J is the 
volume per unit length of fairing, is the x distance from 
the origin to the center of gravity of the fairing, Zis 
the viscous drag force of equation (I-lOe) per unit length of 
fairing segment. is the force on the fairing from the 

cable and p = Since it is more convenient to consider 

the mass of cable and fairing together in the parent problem, 
those terms of equation (I-ll) which will appear directly in 
the inertial and gravitational force terms of the parent 
problem may be lumped in with the cable external force and 
(I-ll) re-written as 

Z = Z '^ + + (Z +(m-pl/)g cos4)*)a + 

c w w a 

(Z +mu*-X*u*)q + (Z*+mx„)q + 2^^^ (1-13) 

q u ^ q G ^ 


Moment equilibrium gives 


-I^q+mXg (w+pv*-qu*) + (m-p I/) gx^ (dcos4)*+iJ)sin(})*) +M^w+M^w+ 


M a+M q+M*q+X*u*^+ZX = 0 
a q^ q^ u p 


(1-14) 


Where is the mass moment of inertia per unit length of the 

fairing about the y axis, and x is the x distance from the 

P 

origin to the point at which the resultant viscous drag force 

acts, a point near the leading edge of the fairing so x will 

P 

be small compared to the chord c. Equation (1-14) may be 
rearranged to give 
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(-Cp^pcXp (a^+a 2 sin<j)*) U^+M^+X^u*)w+ (M^+mx^) w+mx^v*p+ 

(m-pl/) gx-tpsind)* = - (M + (m-p I/) gx^coscp*) a-(M^-mx„u*) q- 
Cj ot vj q u 

(M*-I )q (1-15) 

q y 


Equation (1-15) may now be solved for a as a function 
of ip and w and substituted back into equation (1-13) giving 
finally in terms of ip and w. In view of the fact that the 
solution technique of the parent problem assumes solutions 
sinusoidal in time and that the hydrodynamic derivatives of 
equations (I-lla,b) are functions of frequency, equations (1-13) 
and (1-15) may be solved by assiaming sinusoidal solutions in 
time. Suppose solutions are of the forms 

0 = S 


a = a e 


iwt 


^ — iwt 

ip = ip e 


Z = Z e^“^ 
c c 


Now substituting these solutions into equations (1-13) and 
(1-15) along with the values for the hydrodynamic derivatives 
given by Robinson and Laurmann (1952) leads to 


= [-(l-iH(X) )TTpu*c+(iG(X)-'t)Trpc^ia3-Cpipc(aj^+a2sin<|)*)U^] 
w+[ - (l-^H (A) )iTpu*^c+ (m-pl/)g cos<})*+ ( (-1 +|h (A)+^0 (A ) ) 
Trpu*c^+mu*+iTrpb^u*) ico-( (--^+10 (A) ) 7rpc^+mx^)a)^ ] a 
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[-Cj^^pcXp (aj^+a 2 sin(t)*)U^-T (l-^H(X)) 7Tpu*c^-^TTpb^u* + 

( (.~^+^G (X) ) irpc^+mx^) iw] w+ [mx^v*iw+ (m-pl/) gx^sinct)*] ip = 

(1-17) 

[- (-i (l-^H (X)) iTpu*^c^+ (m-pl/) gx^cosc})*) - ((^ (- 7 +iH (X)) + 
iG(X))TTpu*c^-mx^u*)iu)+ ( (-t47-It (i-^G(X) )iTpc'*-I )a)^] a 

ij y 

where G(X) and H(X) are given in terms of Bessel functions of 
the first and second kinds, that is Jn's and Yn's, as 


H(X) = 2 


J^+Y^+y Jt-J Y, 
o o o 1 o 1 


^ J J,+Y Yt 

G(X) = I- -2-J- ° 1 




(a) C 

where X the reduced frequency is given by X = ^ and the 
argument of the Bessel functions is X/2. We bear in mind that 
the inertial forces lumped into are of 0(a)). 

Equations (1-16) and (1-17) are rather complicated, but 
may be greatly simplified by placing some restrictions on the 
cable fairings and the conditions under which they are towed. 
First assume that the reduced frequency X is such that X<<1. 
This is justified by supposing that the driving frequency for 
lateral motion is the yaw frequency of the towing ship. Then 
typically o)®*! radian per sec. as determined from data presented 
in Korwin-Kroukowsky (1961). Typical fairings usually have a 
chord length less than 1 foot, so for velocities above several 
feet per sec. X<<1. Second, real fairings are usually made 
of materials whose density is close to that of sea water 
so that gravity forces are very small, that is m-pl/=0. With 
little error the fairing geometry may be approximated as 
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shown in figure 1-4. m, 1^, b and may be replaced with 
expressions involving p and c, using the assumption that 
m-pV=0. 

Equations (1-16) and (1-17) may be greatly simplified by 
eliminating inertial, gravitational and thickness effects. 
Assuming low reduced frequency allows us to replace G and H 
by the functions of X to which they tend at low reduced fre¬ 
quency. With the above simplifications a surprising result 
follows. The side force on the cable contributed by the 
hydrodynamics of fairing acting as a thin wing at low reduced 
frequency, completely free to rotate about an axis through 
its leading edge is given by the classical added mass for a 
flat plate of the same span accelerating normal to its plane. 

To this may be added an inertial contribution proportional to 
the rate of change of the kiting angle arising from the center 
of gravity of the fairing not lying at the axis of the 
coordinate system, and the lateral component of viscous drag 
given in equation (I-lOa). 

To elucidate the fairing hydrodynamics the thin wing 
theory may be applied to the case of a massless flat plate to 
which the fairing has been reduced proceeding in the x 
direction with velocity u* being driven in harmonic oscillations 
normal to its line of travel by a force applied at its leading 
edge. The plate is free to rotate about its leading edge 
which is the y axis of a right-handed coordinate system. 
Equations for the forces and moments on the plate are 
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(1-18) 


Z*z + Z-z + Z a + Z*a + Z-ci = f 
z z a a a ex 

M*z + M..Z + M a + M*a + M..a = 0 (1-19) 

z z a a a 

As before z = z a = a e^^^ and f = f e^^^. 

o ' o ex o 

Substituting into (1-18) and (1-19) the appropriate 
expressions for the hydrodynamic derivatives from Robinson 
and Laurmann (1952) gives 

- (l-'^H) Trpu*ciwz^- (^G-^) TTpc^ca^z^- (1-^H) iTpu*^ca^ + 

(-l+'iH+;^G) Trpu*c^ia)a^ + (i-|G) irpc^u^a^ = f^ (1-20) 

-i (1-^H) TTpu*c^ia)Z^+ (i-^G) irpc^w^ z^-^ (1-^H) Trpu*^c^a^ + 
(-•^t'iTH+iG) TTpu*c^ ia)a^-(- tI^+I^G) TTpc‘*a)^a^ = 0 (1-21) 


Equation (1-21) may be solved for giving 

Xz 

-(i(-4+iH) + (i-iG)X]-—^ 

a = -2- (1-22) 

° [ (-^+iH)+i (-l+l-yH+iG) X+ (tIt~'ItG) X^ ] 


where —. = X 

U’^ 

similarly equation (1-20) may be rewritten as 

Xz 

[-i(l-^H)-(iG-t)X]-^ + [-(l-iH) + (-l+^H+iG)iX 
f 


(i-|G)X2]a^ = 


o 

TTPU*'‘^C 


+ 


(1-23) 
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and substituting into (1-23) from (1-22) gives 


(1-24) 

where A(X) = t (-l+^H) -i- (-1 +|h+^G) iX+ (^-|G )\^ ] 

[ (-1+^H)+ (-1+|H+^G) iX+ (^“Ig) ] 

H(X), G(X), and A(X) must now be evaluated for low 

reduced frequency. This is done by looking at series for 

the Bessel functions J , J,, Y , Y, for small values of their 

o 1 o 1 

arguments. From Jahnke and Emde (1945) 

J^(X/2) = 1 - ^ + 0(X‘‘) 

= I + 0(XM 

Yi(X/2) = |[- X I In j + 0(X^ln X)] 

where y = 1.781. Substituting these relationships into the 
expressions for G and H leads to 

H = |2L + |i in2 lA - 3^ ^ O(X^ln^X) 

G = -In ^ + |A In JA + 0(X2ln3X) 
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Substituting these expressions into the expression for A(X) 
gives 


A(X) = 1 - + 0(x2ln2X) 


(1-25) 


Now putting A(X), G(X), and H(X) into equation (1-24) gives 


(| + 0(X2ln2X)) 



f 

o 

Trpu*^c 


or 


f = -^Trpc^ (ioj) ^ z . 

O 4 K X / Q 


(1-26) 


Noticing that io) indicates differentiation with respect to 
time/ the results of equation (1-26) combined with (I-lOa) 
and the term proportional to caused by the center of gravity 
of the fairing not being on the coordinate axis give 

II - If - Cd 

(1-27) 

Then the lateral hydrodynamic force on a segment of fairing 
of length ds is 


Z = -ipCj^cds (aj^+a2sin(t)*)U^ 0 - ^-Trpc^ds 
4mx, 


^ V. as If 


(1-28) 
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Results 


The total hydrodynamic force acting on the cable in the 
X direction may be written by combining equations (I-2a) and 
(I-lOb) giving 

X = -ipirb^ds - Cr^ipcds [ {a,+a„sin(J)*} sin4)*U^ + 

{ (aj^+a2sin4)*) 2sin^4)*+ (aj^+2a2sin4)*) cos^(J)*}U^u 

(1-29) 

+ a, sin(i)*cos4)*U 
1 o 

Similarly the total hydrodynamic force acting on the cable in 
the y direction is given by combining equations (I-2b) and 
(I-9b) 

= -ipiTb^ds u* 1^ - Cj^^pcds [ {bj^+b 2 Cos 4 )* }cos<})*U^ 

+bj^sin(j)*cos4)*U^d+{2 (bj^+b2CoS(t)*) cos^cp* 

+ (bj^+2b2COS(j)*) sin^cj)* }U^^] (1-30) 

and finally the total hydrodynamic force on the cable in the 
lateral direction is given by equation (1-28). 

The above results of the compound hydrodynamic theory of 
the three dimensional unsteady motions of a cable fairing 
segment leave much to be desired in three areas. First the 
theory is more a collection of results bound together by 
assumption than a single cohesive formulation. Second, it 
lacks experimental evidence by which to test its validity. 
Third, the case for the non-neutrally buoyant fairing with 
inertial and thickness effects retained as in equations (1-16) 
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and (1-17) needs to be worked out. The usefulness of the 
theory in the absence of a more comprehensive one is that it 
provides a basis from which to attack the parent problem. 
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Appendix II 

HYDRODYNAMICS OF THE TOWED BODY 

The body under tow in this thesis is a sphere towed from 
its center of gravity and buoyancy. This is, of course, the 
simplest possible towed body and provides the simplest dynamic 
boundary condition at the end of the cable. Real towed bodies 
are often quite complicated,employing fairings to reduce drag, 
depressors to promote deeper towing, and control surfaces to 
alter local orientation, and in very advanced systems, prop¬ 
ulsion. These latter bodies may impose very complicated 
boundary conditions on the cable. 

The equations of motion of the sphere derive from 
considering the inertial forces of its mass and added mass, 
its weight, and its drag all balanced by the tension in the 
cable. Components of the inertial forces in the respective 
coordinate directions are 


F 


(i) 

X 



= -<«B + H' 

= -(Mg + Ug)(|J + u* ||) 


(II-l) 


where is the mass of the body and vi_ is the added mass 

J3 

which for a deeply submerged sphere is 1/2 the displacement. 
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Gravitational forces are 




(Mg - pVg) g(cos(|)* - sin(()*$) 


= -(M„ - pv„) g(sin(f)* + cos(j)*$) (II-2) 

y r> Jd 

= (Mg - pVg) g sin(|)*tl) 

where V„ is the volume of the sphere and p the density of the 

Jd 

water. 

brag forces are found by considering the sphere to be 
moving through the water with velocity W. The drag is then 
given by = -^pC^ S|w|w where is a drag coefficient 

which is a function of Reynold's number and A is the frontal 
area of the sphere. But 


|wl 


/u* ^ + V* ^ 


( 1 + 


u*u v*v . 

u*^+v*^ u* ^+v* ^' 


retaining only the steady terms and the terms linear in the 
perturbation velocities. Drag force components are 


,(d) 

X 

= -hc,^ 


|w| 

(u* 

+ U) 

.(d) 

y 

= -7PC^ 

^B 

|W| 

(v* 

+ V) 

,(d) 

z 

= -rPCrp 

^B 

|w| 

w 



(II-3) 
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And finally 


X 

= ''b ^o 

/U* u 

(ii— + r— + 

o o 

a. 2 ^ 

u* u 

IP" 

o 

U*V*Vx 

o 

>1 

= ^B ^o 

/V* V 

(— + — + 
U 

o o 

v*^v . 

v*u*u^ 

iP"' 

o 

z 

= ^B 

w 



where U 

= /u*'^ + 





o 


Tension acts only in the y direction along the cable and is 
given by 

T* + T (II-5) 


Combining equations (II-l, 2, 4, 5) gives the equations of 
motion for a towed sphere• 

/.. . ^ /9Q 1 ^ tt 2 /U* , u . (u*^u+u*v*v) ^ , 

<Jt " ^ 3t^"^PSVo<U“ 0“ -iP-) 


(Mg-pVg) g(cos(f>* - sin(|)*$) + f^ = 0 


. \ /9V , *9$x 1 ^ , tt 2 /V* , V , (v*^v+v*u*u) ^ 

-'Vb> <3t u~ -iP-> 


-(M„-pV„) g(sin(l>* + cos4)*$) +T* + T+ f = 0 
ID D y 


-(Mg+yg)(|^ + v*||-) +(Mg-pVg) g ip sincp* -ipC^AgU^ w + f^=0 

(II-6) 

In these euqations f , f and f are included as the 
^ X y z 

source of possible perturbing forces such as might be created 
by body control mechanisms. For the simplest case these 
vanish. 
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Considering a general body towed at some arbitrary point 
vice a sphere towed at its center of gravity and buoyancy 
complicates matters. Moments must be considered in writing 
the equations of motion of the body and method of attachment 
of the body to the tow cable becomes a matter of critical 
interest. Eames (1966) considered the general problem in some 
detail. 
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Appendix III 
SPECIAL CASES 

The elements of the coefficient matrix of the cable 
equations of motion are constant for certain special cases. 
Solutions can then be constructed from the eigenvalues and 
eigenvectors of the coefficient matrix. Two cases are con¬ 
sidered. In one case the cable is neutrally buoyant and has 
no viscous drag. The cable tows in a straight line at what¬ 
ever trail angle is specified by the towed body weight and 
drag. As before, the lateral motions may be considered 
separately from the transverse-longitudinal motions. In a 
second case the cable and the towed body are neutrally buoy¬ 
ant. Normal viscous drag acts on the cable but no tangential 
drag. The trail angle is zero, that is the body tows directly 
astern of the tow point. In this case lateral motions, trans¬ 
verse motions and longitudinal motions are all decoupled and 
the lateral and transverse modes are identical in form 
differing only in the value of cable added mass. The point 
is to show that the results obtained from these special cases 
by a simple analysis agree with the results obtained from the 
general formulation and analysis when specialized for the 
same cases. 

As mentioned before in Chapter 3 when the coefficient 
matrix is constant solutions are obtained as follows 

^ = [a] X (III-l) 
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where x is the appropriate vector of n dependent variables 
and [a] is an n x n coefficient matrix. 


|a-Al| = 0 


(III-2) 


gives n eigenvalues of [a]. 



(III-3) 


th th 

gives i^ eigenvector of [a] associated with i eigenvalue. 


^ - \ - r 

X = E c^ e 


i 


(III-4) 


gives the solution in terms of the arbitrary constants of the 


c^. Boundary conditions provide n equations in the n arbitrary 
constants when the equations (III-4) are substituted into them 
for the values of the dependent variables. Solution of these 
equations gives the values of the n arbitrary constants, which 
are substituted back into (III-4) giving the final form of the 
solution. 

The equations of motion are first written in their com¬ 
plex normal forms and specialized for the particular case at 
hand. The above steps are accomplished. The results are 
transformed into the local vertical coordinate system and the 
motions at the towed body are compared with those given by 
the general formulation, method and program. In all cases 
the comparison shows agreement. Even for these simple cases 
the algebra and the computation of the results is tedious so 
simple arithmetic programs were written to do this work. 
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Neutrally Buoyant Cable Without Drag 


Lateral Motions 
The equations of motion are 


rdw] 


cos* 

. ,cos^$ 

N - 

i)' 



iv(i+n) 

N 

ivcos<I>. 

N 

j 



f \ 


w 



j 

J 


(III-5) 


withd=l, D=0, h = 0, 

The steady configuration equations give 

||=0, ^=0, N = C/I+taiP¥, $ = = tan"^ (^) 


(III-6) 


The boundary conditions are 

5 = 0 (iv(M+k)+C)w + (iv (M+k) cos$-(M-a) sin$) = 0 


(111-7) 


C = 1 


w + 


cos$ = e 


The eigenvalues are 


Xi = f ± g 
2 

f = -ivp; g = 



iv/p^+q^ 


/ (1+n) 

/ N6 


Recalling that = w + cos$i|j we may, after some algebra, 
write 
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w = ee 


-ivp (C-1) (“ [ (cos<I)p+^) C+ (M-a) sin<I)p] sin (v/p^+q^^) + 

(- [ (coS(I)p+^) C+ (M-a) sin<I)p] sin (v/p^+q^ ) + 


/p^+q^N^cos$ _ _ 

i[-^- cos (v/p^+q""?) - - - sin (v/p^+^c) ]) 

1 [—*-—-cos(v/p^+q^) - —g- sin(v/p^+q^ )]) 


(III-8) 


from which the magnitude of at ^ = 0 may readily be cal¬ 
culated. For a test cable with M= 3.0, a = 2.0, k = 1.0, 

6 = 0.01, ej^= 1.0, = 0.0, c = 1.0, n = 1.0 agreement 

between the above formulas and the general program was exact 
to at least five significant digits for values of v from 0.0 
to 1.0. 


Transverse-Longitudinal Motions 
The equations of motion are 


fdu' 

dC 


ivcosij) (l+y) 

N 

i rCOS^$ 1 ,x . Nx 1 

IV [ ^ (1+ ) ] 

N 0 Y 

0 

> 

0 


✓ 

U 

d(}) 


iv(1+y) 

ivcos$ 

0 

0 



dC 


N 

N 





dv 


ivsin$(1+y) 

-ivsin'I’cos*!) 

0 

iv 


V 

d? 


N 

y6 




0 

iv (1+y)sin$ 

iv 

0 

J 


n 

1 i 


(III-9) 
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The steady configuration equations are as above for lateral 
motions. The boundary conditions are 

C=0 [iv (M+k)+C (1+sin^ $) ] u+[-iv (M+k) cos'J'+(M-a) sinf] 4i + 
Csin$cos$v = 0 


Csin^cos'J* u + [iv (M+k) sin$+(M-a) cos$] (() + 
[iv(M+k) + C(l+cos^$)]v -n = 0 


(III-IO) 


5=1 u-cos$(() = a = sin$ - cos$ he. 

v+sin#((> = 3 = cos$ + sin$ he 


The eigenvalues are 


. _ ivcos$ . / 2 I 4 (l+y)N(1 +N/y). 

- “In - < ^ -6cos"<I.-^ 


A 3 


±iv 


/ 1 _ 

/ y<5 


(III-ll) 


The solutions are given by 


U 


r > 

ai 


0.2 


' 0 ^ 

; 

' 0 ' 

4> 

_ 

biUi 


CM 

CM 

A 2 C 1 

e 4- 

0 

e^^^+ 

0 

V 


CiUi 


C 2 OL 2 


0 L 3 


a4 

nj 


[dial. 


,d2a2j 


dsOsj 


I^d4 ai, j 




(III-12) 
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where the b's, c's and d's are combinations of eigenvalues and 
coefficient matrix elements found by solving for the eigen¬ 
vectors and the a's are the arbitrary constants of the 
eigenvectors found by substituting equations (III-12) into 
equations (III-IO). The resulting algrbra is tedious so 
results were obtained by computer solutions from this point 
on. The magnitudes of u^, v^ and n at ?=0 were calculated and 
compared with results from the general program for the 
following cable. M=3.0, a=2.0, k=1.0, C=1.0, y= 100.0, 6=0.01, 
y=1.0, 4ap=1.0, -6aj=0.0, h(i^=1.0, hZj=0.0 with v between 0.0 
and 1.0. Results compared favorably through the third 
significant digit. The small remaining difference seems to be 
the result of the matrix inversion carried out in the special 
case solution where solving for the arbitrary constants of the 
eigenvectors by Cramer's Rule. The largest difference occurred 
near a reasonance condition at v=0.3 and was less than 2%. 

Such a reasonance is to be expected since these cases have no 
damping except that which is introduced by the towed body 
boundary condition. 


Neutrally Buoyant Cable and Body With Normal Drag 

Lateral or Transverse Motions 
The equations of motion are 


"dw' 

d? 


(iv(l+n)+D 7 ) 

N 

. ,1 1,1 

(n - 6^ 


< > 

w 

d? 


(iv(l+n)+D 7 ) 

N 

+ 

N 

J 


J 


(III-13) 
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with d=l, N=C, bi=b 2 = 0, M=a and <J>=0 so sin$=0 and cos4>=l and 
Y=100.0. Steady configuration of the cable is trailing dead 
astern. Boundary conditions are 


5 = 0 (“iv (M+k)-C)w - iv(M+k)i/j = 0 


C = 1 


w 


\lJ = 0 


(III-14) 


The eigenvalues are 


Xi = f ± g 
2 


f = 


-D: 

2N 


ivn 

2N 


(III-15) 


g = /dI - v2(ri^+ + i2vD7(ri+ 1^) 


Suppose we also define for convenience three more constants 


_ CD7 _ 2(M+k) 

" 2N “ ^ 6 


B = vc( 25 f + 1) 


(TII-16) 


G = A + iB 


now as in the first case above we may write 


w = ee 
o 


f(5-l) j- GCe^^-e + qC(e^^+e j 

G(e^ -e + gC(e^ +e 


(III-17) 
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The magnitude of at 5=0 may be calculated easily and com¬ 
pared with results for the general programs. Here as before 
the results are exactly similar through the fifth significant 
digit. The equations for the transverse motion are identical 
with n replaced by y and w, ijj replaced by u, (f). The magnitude 
of at 5=0 was calculated and again compared exactly with 
results from the general program through the fifth significant 
digit. 


Longitudinal Motion 
The equations of motion are 


rdv] 

f 

d5 


dn 


dC 

(. j 



0 


iv 


> 

f N 


V 


n 

j 

1 . J 



(III-18) 


for same cable parameters as in the lateral or transverse 
case. The boundary conditions are 


5 = 0 -iv(M+k)v - 2C V + n = 0 


5 = 1 V = B = -iu 


(III-19) 


The eigenvalues are 


X 1 = 
2 



( 111 - 20 ) 
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The solutions are 


V = iu 


[/^ cos 0^^—) - v(M+k)sin( 

/yfi 


/y^ 


) + 2iCsin(^^^)] 
/7T 


[/yT cos ) - V (M+k) sin (—^) + 2iCsin (—^) ] 

/y^ /y6 /y6 

(III-21) 

[/y6[2C cos(^) + i(/yTsin(^)+v(M+k)cos(^))] 
/y^ /yS” 


[/y6 cos (-) - V (M+k) sin (-) +2iC sin (-) ] 

/y6 /yT /y6 

(III-22) 


From these solutions, as before, the magnitudes of u^ 
and n were calculated and compared to the result from the 
general program. Again agreement was exact through five 
significant digits. 

The results of these special cases match closely the 
results of the general program lending credence to the 
niamerical method employed. 
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Appendix IV 
TRIAL CABLE 


The basic towed cable system used in Chapter 4 was 


developed from the following specifications. 

Tow cable: length-1200 ft., diameter-1/2 in., weight 
per foot of cable plus fairing-1/2 lb., 
material-steel with modulus of elasticity 
20 X 10® psi(l) 

Fairing: breadth-0.555 in., mean chord-2.22 in., 

drag coefficient- 0 . 1 (2) 

Body: sphere-40 in. diameter, weight-2482 lbs., 

displacement-1241 lbs., drag coefficient- 
0.15 (3) 


Tow Speed: 16 feet per sec. (approximately 10 kts.) 

Frequency range of excitation: Angular frequency - 

0.0 to 1.0 


The above specifications give the following cable and body 
parameters using the definitions of Chapter 2. 


d=0.125 

b=0.555 


Y=4160.0 
c=2.22 


M=4.14 

y=0.214 


C=0.56 

k=1.03 


a=2.07 D=9.48 

n=3.44 5=0.5 


6=0.0066 

h=0.5 


^^^Meyers, Holm, McAllister, Handbook of Ocean and Underwater 
Engineering , New York, 1969. 

( 2 ) 

Fathom Oceanology Limited Brochure Describing the cable 
fairing called "Flexnose" distributed by Fathom Oceanology 
Limited, 863 Rangeview Road, Port Credit, Ontario. 

(3) 

Hoerner, S.F., Fluid Dynamic Drag , Midland Park, New 
Jersey, 1965. 
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